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, Abstract. We present the X-ray properties and scaling relations of a flux-limited morphology-unbiased sample of 12 X-ray 

^ ■ luminous galaxy clusters at redshift around 0.2 based on XMM-Newton observations. The scaled radial profiles are character- 

' ized by a self-similar behavior at radii outside the cluster cores (> 0.2r5oo) for the temperature (T oc r""^*), surface brightness, 

fT^ , entropy (5 oc r""), gas mass and total mass. The cluster cores contribute up to 70% of the bolometric X-ray luminosity. The 

' X-ray scaling relations and their scatter are sensitive to the presence of the cool cores. Using the X-ray luminosity corrected for 

the cluster central region and the temperature measured excluding the cluster central region, the normalization agrees to better 
^— , than 10% for the cool core clusters and non-cool core clusters, irrelevant to the cluster morphology. No evolution of the X-ray 

scaling relations was observed comparing this sample to the nearby and more distant samples. With the current observations, 
the cluster temperature and luminosity can be used as reliable mass indicators with the mass scatter within 20%. Mass discrep- 
. ancies remain between X-ray and lensing and lead to larger scatter in the scaling relations using the lensing masses (e.g. ~ 40% 

CL(' for the luminosity-mass relation) than using the X-ray masses (< 20%) due to the possible reasons discussed. 

2 ' Key words. Cosmology: observations - Galaxies: clusters: general - X-rays: galaxies: clusters - (Cosmology:) dark matter 
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1. Introduction 2006) and galaxy clusters combining the Sunyaev-Zeldovich 

effect (Sunyaev & Zeldovich 1972, the SZ effect) and ther- 



Both the gravitational growth of density fluctuations and the 



mal Bremsstrahlung X-ray approaches (Molnar et al. 2002), in 



$H ' expansion history of the Universe can be used to constrain the j i j j 

rri • , which the luminosity distance and volume are used to deter- 

>r cosmology. ■ »u ■ i 

mine the metric, respectively. 

The gravitational growth of the fluctuation can be mea- 
sured by, e.g. the evolution of the galaxy cluster mass func- The multi-pole structure of the cosmic microwave back- 
tion (e.g. Schuecker et al. 2003). The most massive clusters ground (CMB) anisotropy power spectrum depends on the nor- 
show the cleanest results in comparing theory with observa- malized overall amount of dark matter (DM) in the Universe, 
tions. The mass function of luminous galaxy clusters probes The WMAP three year results imply values of the parame- 
the cosmic evolution of large-scale structure (LSS) and is thus ters Qm/;^, Qb/i^, h, n^, and erg of 0.127+° °°3, 0.0223;^Q °Q°g, 
an extremely effective test of cosmological models. It is sen- 0.73+q;o3, 0.95 1+°;°}^, and 0.74+° °^ (Spergel et al. 2006). 
sitive to the matter density, Q.^, and the amphtude of the cos- Combining the CMB approach with the other cosmological 
mic power spectrum on cluster scales, erg (e.g. Schuecker et al. tests, e.g. galaxy cluster surveys, the large-scale galaxy distri- 
2003). bution and/or SN la, is sensitive to the effects of dark energy, 

The expansion history of the Universe can be used to deter- the density of which is characterized by the parameter Qa and 

mine the metric and thus to constrain the cosmological parame- its time evolution by the equation of state parameter w{z) (e.g. 

ters. Typical examples are Supernova la (SN la, e.g. Astier et al. Vikhlinin et al. 2002; Allen et al. 2004; Chen & Ratra 2004; 

Hu et al. 2006). 

Send offprint requests to: Y.-Y. Zhang 

* This work is based on observations made with the XMM-Newton, The construction of the mass function of galaxy clusters for 

an ESA science mission with instruments and contributions directly large cosmological cluster samples is based on the cluster X- 

funded by ESA member states and the USA (NASA). ray temperature/luminosity estimates and the mass-observable 
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scaling relations (Reiprich & Bohringer 2002). A better un- 
derstanding of the mass-observable scaling relations and their 
scatter is of prime importance for galaxy clusters as a unique 
means to study the LSS and thus to constrain the cosmologi- 
cal parameters (Voit 2005). Precise mass estimates and mass- 
observable scaling relations can potentially push the current 
measurement uncertainty (10-30%) of the cosmological tests 
to higher precision (Smith et al. 2003, 2005; Bardeau et al. 
2006). 

The mass distributions in galaxy clusters can be probed in a 
variety of ways using: (i) the X-ray gas density and temperature 
distributions under the assumptions of hydrostatic equilibrium 
and spherical symmetry, (ii) the velocity dispersion and spa- 
tial distribution of cluster galaxies, (iii) the distortion of back- 
ground galaxies caused by the gravitational lensing effect, and 
(iv) the SZ effect. For the relaxed low-redshift clusters, the un- 
certainties are small in the masses measured from the precise 
intra-cluster medium (ICM) temperature and density profiles 
up to large radii (e.g. rsoo) using X-ray observations of XMM- 
Newton and Chandra (Allen et al. 2004; Pointecouteau et al. 
2005; Arnaud et al. 2005; Vikhlinin et al. 2006a). Particularly, 
the X-ray results for the relaxed clusters have achieved a very 
high level of precision (< 10%) in Arnaud et al. (2005) and 
Vikhlinin et al. (2006a). The methods to measure cluster mass 
using the velocity dispersion of cluster galaxies and gravita- 
tional lensing are improving (Girardi et al. 1998; Smith et al. 
2003, 2005; Bardeau et al. 2006). Up to now, it is technically 
more difficult for the SZ effect to resolve radial temperature 
distributions in galaxy clusters to obtain precise cluster masses 
but the future holds promise (Zhang & Wu 2003). 

Some interpretations of the cluster dynamics, such as spher- 
ical symmetry and hydrostatic equilibrium, are required to 
obtain the cluster masses for the X-ray approach (also for 
the SZ effect). Substructure is observed in galaxy clusters 
and the frequency of its occurrence in ROSAT observations 
has for example been estimated to be of the order of about 
52 + 7% (Schuecker et al. 2001). Using high quaUty XMM- 
Newton data, Finoguenov et al. (2005) found that substructures 
(on > 10% level) can be observed in all the clusters in the 
REFLEX-DXL sample at z ~ 0.3. Based on the strong lens- 
ing results using high resolution HST data. Smith et al. (2005) 
measured the substructure fractions and identified the clusters 
with multi-modal DM morphologies. The up-coming XMM- 
Newton/Chandra results of comprehensive samples without ex- 
clusion of the mergers (HIFLUGS, Reiprich & Bohringer 2002, 
Reiprich 2006; LP, Bohringer et al. 2007), would be promising 
to investigate the mass measurement accuracy and the scaling 
relations for a morphology-unbiased sample at low redshifts. 

Independent of the X-ray method, gravitational lensing pro- 
vides a direct measure of projected cluster masses irrespective 
to their dynamical state. Comparison of ground-based lens- 
ing data including estimated arc redshifts with ASCA/ROSAT 
X-ray data show that the masses between X-ray and lens- 
ing are closely consistent for relaxed clusters (Allen 1998). 
However, detailed investigation with high quality (HST/XMM- 
Newton/Chandra) data reveals that for a few apparently "re- 
laxed" clusters in X-rays there remain large mass discrepan- 
cies between X-ray and lensing, especially on small scales, e.g. 



CL0024+17 (Kneib et al. 2003; Ota et al. 2004; Zhang et al. 
2005a) and Abell 2218 (Soucail et al. 2004; Smith et al. 2005; 
Pratt et al. 2005). Stanek et al. (2006) illustrate the interplay 
between parameters and sources of systematic eiTor in cosmo- 
logical applications using galaxy clusters and stress the need 
for independent calibration of the L-M relation such as grav- 
itational weak lensing and X-ray approaches. It is extremely 
important to investigate cluster dynamics and to calibrate the 
scaling relations of a set of morphology-unbiased cluster sam- 
ples at different redshifts to improve the understanding of in- 
trinsic scatter and to control the sources of systematic eiTors. 

Recently two morphology-unbiased flux-limited samples at 
medium distant redshifts were constructed. One (Smith et al. 
2005; Bardeau et al. 2005, 2006; hereafter the pilot LoCuSS 
sample) is in the redshift range 0.176 < z < 0.253 with 
L > 8.0 X 10'*^ erg s ' based on the XBACS in EbeHng et al. 
(1996). The other (Zhang et al. 2004a, 2006; Finoguenov 2005; 
Bohringer et al. 2006; the REFLEX-DXL sample) is in the red- 
shift range 0.258 < z < 0.308 with L > 5.9xl0'*4ergs-i based 
on the REFLEX survey in Bohringer et al. (2004). The selec- 
tion effect is well understood for both samples (Smith et al. 
2005; Bohringer et al. 2006), consisting of massive galaxy clus- 
ters spanning the full range of cluster morphologies (Bardeau 
et al. 2006; Zhang et al. 2006). Using the same method de- 
scribed in Zhang et al. (2006) for the REFLEX-DXL sam- 
ple, we present here the studies of the pilot LoCuSS sam- 
ple of 12 clusters (Abell 2219 is flared) based on high qual- 
ity XMM-Newton observation. Smith et al. (2005) performed 
the strong lensing studies for most clusters in this sample ex- 
cept for Abell 1689 (Broadhurst et al. 2005a; Halkola et al. 
2006) and Abell 2667 (Covone et al. 2006). The strong lensing 
studies of Abell 2218 can also be found in Kneib et al. (1996, 
2004). Weak lensing analysis using HST, CFH12k and Subaru 
can be found in Bardeau et al. (2005, 2006) and Broadhurst et 
al. (2005b, for Abell 1689). We note that Abell 773 has no weak 
lensing mass estimates because no CFH12k data were taken for 
this cluster. We will compare the X-ray masses with the lensing 
masses (e.g. Kneib et al. 1996, 2004; Smith et al. 2003, 2005; 
Bardeau et al. 2006) to discuss possible mass discrepancies. 

The main goals of this work are: ( 1 ) to derive precise cluster 
mass and gas mass fraction, (2) to investigate the self-similarity 
of the scaled profiles of the X-ray properties, (3) to seek for 
proper definition of the global quantities for the X-ray scal- 
ing relations with the least scatter, (4) to study the evolution of 
the X-ray scaling relations combining the pilot LoCuSS sample 
with the nearby and more distant samples, and (5) to compare 
the X-ray and lensing masses and to understand the scatter in 
the luminosity-mass scaling relations using the lensing masses. 

The data reduction is described in Sect. |2] and the X-ray 
properties of the ICM of the sample in Sect. |3] We investigate 
the self-similarity of the scaled profiles of the X-ray properties 
in Sect, m and the X-ray scaling relations in Sect.|5] In Sect.|6] 
we compare the X-ray and lensing masses, discuss the pecu- 
liarities in individual clusters, and address the possible expla- 
nation for the scatter in the relation between the lensing mass 
and X-ray luminosity. We our the conclusions in Sect.|7] Unless 
explicitly stated otherwise, we adopt a flat ACDM cosmology 
with the density parameter = 0.3 and the Hubble con- 
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slant Hq - 70 km s"' Mpc"^ We adopt the solar abundance 
table of Anders & Grevesse (1989). Confidence intervals cor- 
respond to the 68% confidence level. We apply the Orthogonal 
Distance Regression method (ODRPACK 2.010, e.g. Boggs et 
al. 1987) taking into account measurement errors on both vari- 
ables to determine the parameters and their error bars of the 
fitting throughout this paper. We use Monte Carlo simulations 
for the uncertainty propagation on all quantities of interest. 

2. Data reduction 

2.1. Data preparation 

All 13 clusters of the pilot LoCuSS sample at z ~ 0.2 were 
observed by XMM-Newton. However, the XMM-Newton ob- 
servations of Abell 2219 are seriously contaminated by soft 
proton flares and are therefore excluded in this work. In total, 
12 galaxy clusters observed by XMM-Newton were uniformly 
analyzed using the same method as described in Zhang et al. 
(2006). We use the XMMSAS v6.5.0 software for the data re- 
duction. In Table|AT|(5ee the electronic edition of the Journal) 
we present a detailed XMM-Newton log of the observations for 
all clusters. 

All observations were performed with medium or thin filter 
for three detectors. The MOS data were taken in Full Frame 
(FF) mode except for MOSl of Abell 1835. The pn data were 
taken in Extended Full Frame (EFF) mode or FF mode. For 
pn, the fractions of the out-of-time (OOT) effect are 2.32% and 
6.30% for the EFF mode and FF mode, respectively. An OOT 
event file is created and used to statistically remove the OOT 
effect. 

Above 10 keV (12 keV), there is little X-ray emission from 
clusters detected with MOS (pn) due to the low telescope ef- 
ficiency at these energies. The particle background therefore 
dominates. The light curve in the energy range 10-12 keV (12- 
14 keV) for MOS (pn), binned in 100 s intervals, is used to 
monitor the particle background and to excise periods of high 
particle flux. Since episodes of "soft proton flares" (De Luca & 
Molendi 2004) were detected in the soft band, the light curve 
of the 0.3-10 keV band, binned in 10 s intervals, is used to 
monitor and to excise the soft proton flares. A 10 s interval 
bin size is chosen for the soft band to ensure a similar good 
photon statistic as for the hard band. The average and variance 
of the count rate have been interactively determined for each 
light curve from the count rate histogram. Good time intervals 
(GTIs) are those intervals with count rates below the thresh- 
old, which is defined as 2-cr above the quiet average. The GTIs 
of both the hard band and the soft band are used to screen the 
data. The background observations are screened by the GTIs 
of the background data, which are produced using exactly the 
same thresholds as for the corresponding target observations. 
Settings of FLAG = and PATTERN < 13 (PATTERN < 5) 
for MOS (pn) are used in the screening process. 

An "edetect_chain" command was used to detect point-like 
sources. Point sources were subtracted before the further data 
reduction. 

' http://www.netlib.org/odrpack/ and references here 



2.2. Background subtraction 

The background consists of several components exhibiting dif- 
ferent spectral and temporal characteristics (e.g. De Luca & 
Molendi 2001; Lumb et al. 2002; Read & Ponman 2003). The 
background components can approximately be divided into two 
groups (e.g. Zhang et al. 2004a). One contains the background 
components showing significant vignetting, e.g. the cosmic X- 
ray background (CXB). The other contains the components 
with very little or no vignetting, e.g. particle-induced back- 
ground. 

Suitable background observations guarantee similar back- 
ground components showing vignetting as for the target in the 
same detector coordinates. We choose the blank sky accumula- 
tions in the Chandra Deep Field South (CDFS) as background. 
The background observations were processed in the same way 
as the cluster observations. The CDFS observations used the 
thin filter for all detectors and the FF/EFF mode for MOS/pn. 
The medium filter was used for some target observations. Using 
the medium filter the background is different from the back- 
ground using the thin filter at energies below 0.7 keV. Therefore 
we performed all the analysis at energies above 0.7 keV, in 
which the difference of the background is negligible. One can 
subtract the background extracted in the same detector coor- 
dinates as for the target. The cluster emission covers the in- 
ner part of the field of view (FOV), R < 8', and leaves the 
outer region to exam the background. The difference between 
the background in the target and background observations is 
taken into account as a residual background in the background 
subtraction. We applied such a double-background subtraction 
method specially developed for the medium distant clusters in 
Zhang et al. (2006) for spectral imaging analysis. More details 
about this method can be found in Zhang et al. (2006). Note 
independently a similar type of double background subtrac- 
tion was earlier described in Pratt et al. (2001) and formalised 
in Arnaud et al. (2002). Here we only briefly described our 
double-background subtraction method (Zhang et al. 2006) as 
follows. 

2.2.1. Spectral analysis 

The spectra are extracted using the annuli within the trunca- 
tion radius (see rt in Sect. 13.31 and Table [T]) corresponding to 
a S'/ZV of 3 of the observational surface brightness profile. The 
width of each annulus is selected to achieve precise temper- 
ature estimates. In order to obtain temperature measurements 
with uncertainties of ~ 15%, we used the criterion of ~ 2000 
net counts per bin in the 2-7 keV banc0. 

For a given target spectrum (5 a) extracted from the annu- 
lus region with an area of A^ in the target observations, the 
background spectrum is extracted in the background obser- 
vations in the same detector coordinates as for the target spec- 
trum. Hereafter we call the count rate ratio of the target and 
background observations limited to 10-12 keV (12-14 keV) 
for MOS (pn) as cab- The first-order background spectrum 

- The 2-7 keV band is sensitive to the determination of the tem- 
peratures for massive clusters which have cluster temperatures higher 
than 3 keV. 



4 



Zhang et al.: Scaling Relations and Mass Calibration of the X-ray Luminous Galaxy Clusters at z=0.2 



is CAB^a- Because the cluster emission covers radii smaller 
than 8' as shown in Sect. 13.31 (also see the truncation radius 
in Table [1]), the second-order background spectrum can be 
prepared as follows. A spectrum (^o) is extracted from the 
outer region (e.g. 8' < /? < 8.33') with an area of in the 
target observations and its background spectrum (Bo) in the 
same detector coordinates but in the background observations. 
The second-order background spectrum is ^{So - cabBq)- 
Assuming that the second-order background spectrum consists 
of the vignetted components such as CXB, the cluster spec- 
trum is then given by - cp^^^B^ - ^{So - cabBq)- We ap- 
plied the second-order background spectrum with and with- 
out vignetting coiTection in the spectral analysis, and found 
that the difference in temperature and metallicity measure- 
ments increases with radius only up to 3% and 17%, respec- 
tively. Therefore the influence on the spectral measurements 
due to possible non-vignetted components such as the particle- 
induced background in the second-order background is negli- 
gible. 

We performed spectral fitting of the MOS and pn data si- 
multaneously with the XSPEC v 11. 3.1 software. Both the re- 
sponse matrix file (rmf) and auxiliary response file (arf) are 
used to recover the correct spectral shape and normalization 
of the cluster emission component. The following are usually 
taken into account for the rmf and arf, (i) a pure redistribution 
matrix giving the probability that a photon of energy E, once 
detected, will be measured in data channel PI, (ii) the quantum 
efliciency (without any filter, which, in XMM-Newton calibra- 
tion, is called closed filter position) of the CCD detector, (iii) 
filter transmission, and (iv) geometric factors such as bad pixel 
corrections and gap corrections (e.g. around 4% for MOS). The 
vignetting coiTection to effective area for ofF-axis observations 
can be accounted for in the event lists by a weight column cre- 
ated by "evigweight". All spectra are extracted considering the 
vignetting correction by the weight column in the event list pro- 
duced by "evigweight". The on-axis rmf and arf are co-created 
to account for (i) to (iv). We fixed the redshift to the peak value 
of the cluster galaxy redshift histogram (e.g. Smith et al. 2005; 
Bardeau et al. 2006; Covone et al. 2006) and the Galactic ab- 
sorption to the weighted value in Dickey & Lockman (1990). A 
combined model "wabs*mekal" is then used with the on-axis 
arf and rmf in XSPEC for flie fitting. 

2.2.2. Image analysis 

The 0.7-2 keV band is used to derive the surface brightness 
profiles. This ensures an almost temperature-independent X- 
ray emission coefficient over the expected temperature range. 
The vignetting correction to effective area is accounted for by 
the weight column in the event lists created by "evigweight". 
Geometric factors such as bad pixel corrections are accounted 
for in the exposure maps. The width of the radial bins is 2". An 
azimuthally averaged surface brightness profile of the CDFS is 
derived in the same detector coordinates as for the target. The 
count rate ratios of the target and CDFS in the 10-12 keV band 
and 12-14 keV band for MOS and pn, respectively, are used 
to scale the CDFS surface brightness. The residual background 



in each annulus of the surface brightness is the count rate in 
the 0.7-2 keV band of the area scaled residual spectrum ob- 
tained in the spectral analysis. Both the scaled CDFS surface 
brightness profile and the residual background are subtracted 
from the target surface brightness profile. The background sub- 
tracted and vignetting corrected surface brightness profiles for 
three detectors are added into a single profile, and re-binned 
to reach a significance level of at least 3-cr in each annulus. 
The truncation radii rt are listed in Table [T] For most clusters, 
the particle-induced background varies by less than 10% com- 
paring the background observations to the target observations. 
Therefore the dispersion of the re-normalization of the back- 
ground observations is typically 10%. We take into account a 
10% uncertainty of the scaled CDFS background and residual 
background. 

2.3. PSF and de-projection 

In the imaging analysis, we coiTect the PSF effect by fitting the 
observational surface brightness profile with a surface bright- 
ness profile model convolved with the empirical PSF matrices 
(Ghizzardi 2001). 

Using the XMM-Newton point-spread function (PSF) cal- 
ibrations in Ghizzardi (2001) we estimated the redistribution 
fraction of the flux. We found 20% for bins with width about 
0.5' and less than 10% for bins with width greater than 1' ne- 
glecting energy dependent effects. We thus require an annulus 
width greater than 0.5' in the spectral analysis. For such distant 
massive clusters, the PSF effect is only important within 0.5', 
which corresponds < O.lrgoo and introduces an uncertainty to 
the final results of the temperature profiles. This has to be inves- 
tigated using deeper exposures with better photon statistic. The 
PSF blurring can not be completely considered for the spectral 
analysis the same way as for the image analysis because of the 
limited photon statistic. 

The projected temperature is the observed temperature 
from a particular annulus, containing in projection the fore- 
ground and background temperature measurements. Under the 
assumption of spherical symmetry, the gas temperature in each 
spherical shell is derived by de-projecting the projected spectra. 
In this procedure, the inner shells contribute no flux to the outer 
annuli. The projected spectrum in the outermost annulus is thus 
equal to the spectrum in the outermost shell. The projected 
spectrum in the neighboring inner annulus has contributions 
from all the spectra in the shell at the radius of this annulus and 
in the outer shells (e.g. Suto et al. 1998). We de-projected the 
spectra as described in Jia et al. (2004) and performed the spec- 
tral fitting of the de-projected spectra in XSPEC to measure the 
radial temperature and metallicity profiles. 



3. X-ray properties 

The primary parameters of all 13 galaxy clusters are given in 
Table □ 
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3.1. Density contrast 

To determine the global cluster parameters, one needs a fiducial 
outer radius that was defined as follows. The mean cluster den- 
sity contrast is the average density with respect to the critical 
density, 

3M(< r) 



A(< r) 



(1) 



The critical density at redshift z is pdz) - pcoE^(z), where 
E^iz) = Qm(l + z)^ -H Qa -H (1 - Q„ - Qa)(1 + z)2. ta is the 
radius within which the density contrast is A. Ma is the total 
mass within r^. For A = 500, rsoo is the radius within which 
the density contrast is 500 and Mgoo is the total mass within 
r5oo- For A = 2500, r25oo is the radius within which the density 
contrast is 2500 and M2500 is the total mass within r250(). 

3.2. Temperature profiles 

Temperature profiles probe the thermodynamical history of 
galaxy clusters. XMM-Newton (also Chandra), in contrast to 
earlier telescopes, has a less energy-dependent, smaller PSF, 
more reliable to study cluster temperature profiles. We de- 
projected the spectra as described in Sect. 12.31 and performed 
the spectral fitting in XSPEC to obtain the radial measurements 
of the temperature and metallicity. In the left panels of Fig. [T] 
(Also see Fiss. \B.l\\B.5\ in the electronic edition of the Journal 
for the whole sample) we show the radial temperature profiles. 
The temperature profiles are approximated by 

T{r) = exp[-(r - T,flT2\ + Te{\ + r^lTf)-^^ + Tj , (2) 

where T,-, / = 1, ...,7, are simply for parameterization without 
physical meaning. 

A cool core cluster (CCC) often shows both a peaked sur- 
face brightness profile and a steep temperature drop towards the 
cluster center. The classical mass deposition rate is derived by 
the radiative cooling model linking to the X-ray luminosity and 
thus the surface brightness profile (e.g. Peres et al. 1998). Chen 
et al. (2007) defined CCCs by two criteria, (1) scaled mass 
deposition rate M/M500 > lO^'-'/yr, and (2) mass deposition 
rate M > O.OlMo/yr. They found that 60% of the HIFLUGS 
sample are CCCs. With high resolution Chandra observations 
Vikhlinin et al. (2006b) defined CCCs by the cuspness of the 



surface brightness profile (a 



t/logPg 
J log r 



> 0.5 at r = 0.04/;50o), 



and found that 2/3 of the clusters at z < 0.5 and 15% at z > 0.5 
are CCCs, respectively. At z > 0.5, they found no pronounced 
CCCs with a > 0.7. 

For Abell 383, Abell 1835, Abell 2390 and Abell 2667, the 
temperatures drop to the values less than 70% of the maximum 
temperatures towards the cluster center. There is no strict divi- 
sion between CCCs and non-CCCs. We empirically define the 
CCC as the cluster whose temperature drop is greater than 30% 
of the peak temperature towards the cluster center. Hereafter 
these 4 clusters are called CCCs in the sample. As we show 
later, the 4 CCCs show shorter central cooling time, larger cool- 
ing radii (scaled to r5oo) and lower central entropies compared 
to the non-CCCs. Using our definition of the CCC, the frac- 
tions of the CCCs are 33% and 15%, respectively, for the sam- 
ples at z ~ 0.2 in this work and z ~ 0.3 in Zhang et al. (2006). 



This could indicate the evolution of the cool cores for massive 
clusters between z ~ 0.2 and z ~ 0.3 that the CCC fraction 
decreases with increasing redshift. This evolution could have 
already become important at z ~ 0.3. 

However, the Poisson noise is quite significant for such 
small samples. Using re-randomization of the Monte Carlo 
simulated temperature profiles, the CCC fraction is 25^0% for 
the pilot LoCuSS sample, and 15-30% for the REFLEX-DXL 
sample, respectively. The difference in the fraction of the CCCs 
is actually not significant between the samples at redshift ~ 0.2 
and 0.3. The evolution of the cool cores is thus not justified by 
the low significance of our result. Furthermore, the bootstrap 
effect accounting for missing clusters due to (1) the selection 
effect of the sample, and (2) the flares the XMM-Newton ob- 
servational run, can only be recovered by mock catalogs in sim- 
ulations. Additionally, the detection of the temperature drop in 
the cluster center becomes difficult for distant clusters, which 
can enhance the phenomenon that the CCC fraction decreases 
with the increasing redshift. A less resolution-dependent way 
to investigate the evolution of the cool cores is to use the lu- 
minosity based on the flux images, which will be addressed in 
SectlO 

3.3. Surface brightness 

A p model (e.g. Cavaliere & Fusco-Femiano 1976; Jones & 
Forman 1984) is often used to describe electron density pro- 
files in clusters. To obtain an acceptable fit for all clusters in this 
sample, we adopt a double-y6 model of the electron number den- 
sity profile, ne(r) = neOi(l + rVr^j)"'^^' + "002(1 + r^lrl^Y^^I^, 
where n^Q - Meoi + ne02 is the central electron number den- 
sity (see Table |2|i, /? the slope parameter, and r^i and r^2 the 
core radii of the inner and outer components, respectively (e.g. 
Bonamente et al. 2006). The soft band (e.g. 0.7-2 keV) X-ray 
surface brightness profile model S (R), in which R is the pro- 
jected radius, is linked to the radial profile of the ICM electron 
number density n^-ir) as an integral performed along the line- 
of-sight for hot clusters [T > 2 keV), 



Sx(R) ' 



Jr 



d€. 



(3) 



We fitted the observed surface brightness profile by this integral 
convolved with the PSF matrices (middle panels in Fig.lT] also 
see Figs. \B.1HB.5\ in the electronic edition of the Journal for 
the whole sample) and obtained the parameters of the double- 
ts model of the electron density profile. The fit was performed 
within the truncation radius (rt, see Table [TJ corresponding to 
a 5 /A^ of 3 of the observational surface brightness profile. The 
truncation radii r^ are above rsoo for all 12 clusters. The cluster 
cores, referring to the central /? model of the double-y6 model, 
span a broad range up to 0.2 rsot). The cluster cores for the 
CCCs are resolved with the current observations. 

3.4. Mass distribution 

We assume that the ICM is in hydrostatic equilibrium within 
the gravitational potential dominated by DM and its distribu- 
tion appears spherical symmetry. The ICM can then be used to 
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trace the cluster potential. The cluster mass is calculated from 
the X-ray measured ICM density and temperature distributions, 

1 d[n,{r)T{r)] ^ GM« r) 

where // = 0.62 is the mean molecular weight per hydrogen 
atom. Following the method in Zhang et al. (2006), we use a set 
of input parameters of the approximation functions, in which 
jS, MeOi, '"ci (' = 1,2) represent the electron number density pro- 
file ne(r) and Ti (i = 1, ...,7) represent the temperature profile 
T(r), respectively, to compute the mean cluster mass. The mass 
uncertainties are propagated using the uncertainties of the elec- 
tron number density and temperature measurements by Monte 
Carlo simulations as follows. For each cluster, we simulated 
electron density and temperature profiles of a sample of 100 
clusters using the observed electron density and temperature 
profiles and their errors. The mass profiles and other properties 
of the 100 simulated clusters were calculated to derive the error 
bars. 

The observed mass profile was used to estimate M500 and 
rsoo (see Table |2]). The NFW model (e.g. NavaiTO et al. 1997, 
NFW) does not provide an acceptable fit for the observed mass 
profiles. We adopt the acceptable fit of a generalized NFW 
model (e.g. Hernquist 1990; Zhao 1996; Moore et al. 1999, 
Navarro et al. 2004), pouir) = Ps('-/'-s)""'(l + r/r,)°'-\ where 
Ps and rs are the the characteristic density and scale of the halo, 
respectively (see Table O. Note a - Q could be an artificial 
fact in fitting the central data points of the mass profile due to 
low spatial resolution and/or significant PSF effect. In the right 
panels of Fig.[T](A/50 see Fiss. \B.l\\B.5\ in the electronic edition 
of the Journal for the whole sample), we present the observed 
mass profiles and their best generalized NFW model fits. We 
note that the data points in the mass profiles are not completely 
independent, which can introduce uncertainties in the fitting. 

3.5. Gas mass fraction distribution 

The gas mass fraction is an important parameter for cluster 
physics, e.g. heating and cooling processes, and cosmological 
applications using galaxy clusters (e.g. Vikhlinin et al. 2002; 
Allen et al. 2004). The gas mass fraction distribution is defined 
as /gas(< '") = A/gas(< r)/M(< r). The gas mass, total mass and 
gas mass fraction at rsoo are given in Table |2] We obtained an 
average gas mass fraction of 0.13 + 0.01 at rsoo for the pilot 
LoCuSS sample. This result agrees with the average gas mass 
fraction of 0.1 16 + 0.007 at rjoo for the REFLEX-DXL sample 
at z ~ 0.3 in Zhang et al. (2006). The average gas mass frac- 
tion at r25(K) of the pilot LoCuSS sample is 0.107 + 0.005, in 
good agreement with the value of the REFLEX-DXL sample 
in Zhang et al. (2006) giving 0.100 + 0.004, and the value in 
Allen et al. (2002) based on Chandra observations of 7 clusters 
giving/g,s~0.105-0.138/i;o'^'. 

As expected, the gas mass fraction distributions of all the 
clusters are lower than the universal baryon fraction, /b = 
Qb/Qm = 0.176 + 0.019, based on the WMAP three year re- 
sults in Spergel et al. (2006). This is because that the baryons 
in galaxy clusters reside mostly in hot gas together with a small 
fraction of stars as implied in simulations (15% /gas, e.g. Eke et 



al. 1998, Ki-avtsov et al. 2005; 5-7% /gas, e.g. Eke et al. 2005). 
In principle, Qm can be determined from the baryon fraction, 
/b = /eas + /eai, in which a contribution from stars in galaxies 
is given by /gai = 0.02 + O.Ol/ig,,' (White et al. 1993). The gas 
mass fractions of the pilot LoCuSS sample support a low mat- 
ter density Universe as also shown in recent studies (e.g. Allen 
et al. 2002; Ettori et al. 2003; Vikhlinin et al. 2003). 

3.6. Global temperature and metallicity 

Vikhlinin et al. (2005) used the volume average of the radial 
temperature profile in a certain radial range as the global tem- 
perature. We seek well defined global temperatures for the pi- 
lot LoCuSS clusters as follows. We calculated the volume av- 
erage of the radial temperature profile in the radial range of 
O.lrsoo < r < 0.5r50() and 0.2r5oo < r < 0.5r5oo, respectively. 
We also measured the spectral temperatures (T'^p'^'^) and metal- 
licities (Z) in the annuli ofR< 2/3rt, O.lrsoo < R < 0.5r5oo and 
0.2r5oo < R < 0.5r5()o, respectively (see Table |4|i. The spec- 
tral temperature agrees with the volume average of the radial 
temperature profile better in the annulus of 0.2-0.5r50() than in 
the annulus of 0.1-0.5r5oo because the measurements limited 
to 0.2-0.5r5oo are less affected by the cool cores for the CCCs. 

We found an average of 0.29 + 0.02 Zq for the global metal- 
licities of the sample measured in the < 2/3rt regions. For the 8 
non-CCCs of this sample the average metallicity, 0.26+0.02Zq, 
agrees with the value of 0.24 + 0.03 Zq for the REFLEX- 
DXL sample (Zhang et al. 2006), and O.21+[J i°Z0 for 18 distant 
(0.3 < z < 1.3) clusters in Tozzi et al. (2003) within the mea- 
surement uncertainties. The average metallicity of the CCCs of 
the sample, O.36+O.O3Z0, agrees with the value of 0.38+0.07Zo 
for 21 CCCs in Allen & Fabian (1998). 

4. Self-similarity of the scaled profiles of the X-ray 
properties 

Simulations (e.g. NavaiTO et al. 1997, 2004) suggest a self- 
similar structure for galaxy clusters in hierarchical structure 
formation scenarios. The scaled profiles of the X-ray proper- 
ties and their scatter can be used to quantify the structural vari- 
ations. This is a probe to test the regularity of galaxy clusters 
and to understand their formation and evolution. The accuracy 
of the determination of the scaling relations, limited by how 
precise the cluster mass and other global quantities can be es- 
timated, is of prime importance for the cosmological applica- 
tions of clusters of galaxies. 

Because the observational truncation radii (rt) in the surface 
brightness profiles are above rgoo for all clusters but below r2oo 
for most clusters in the sample, we use rgoo for radial scaling. 

The following redshift evolution corrections (e.g. Zhang et 
al. 2006) are usually used to account for the dependence on the 
evolution of the cosmological parameters, 

Sx-E-'(z) (A,,/A,,o)-'-''^/(r), 

S ■ E''\z) (A,,,/A,,o)2/3 <x /(D, 

L-E-'(z) (A,,,/A,,o)-°-5 (X /(D, 

M • E(z) (AJAcfif' oc fiT), 

Mgas • E{z) (A,,JA,,o)0-5 oc f(T), 

r ■ E{z) (A,,,/A,,o)0-5 <x f(T), 
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where Ac,z = IS/r^ + 82(Qni,z - 1) - 39(Qm,z - l)Ms the 
analytic approximation derived from the top-hat spherical col- 
lapse model for a flat Universe (Bryan & Norman 1998) and 
Qni,z the cosmic density parameter at redshift z- The function 
f{T) denotes the best fitting power law parameterization. 

We investigated the self-similarity of the scaled profiles of 
the X-ray properties for this sample as follows. 



4.1. Scaled temperature profiles 

We scaled the radial temperature profiles by ro.2-o.5r5oo ™d rsoo 
(Fig. 12]). Within 0.2r5oo, we observed a temperature drop to at 
least 70% of the maximum value towards the cluster center for 
4 clusters (defined as the CCCs) and an almost constant tem- 
perature distribution for the non-CCCs. Previous observations 
have shown that temperature measurements on scales below 
0.2r5()() tend to show peculiarities linked to the cluster dynam- 
ical history. For example, the temperatures of merger clusters 
can be boosted (Smith et al. 2005). However, the boosting is 
serious mainly in the cluster cores. Using global temperatures 
excluding the < 0.2r5oo regions can thus (1) tight the scaled 
profiles of the X-ray properties, (2) minimize the scatter in the 
X-ray scaling relations and (3) reach a better agreement in the 
normalization between the X-ray scaling relations for the CCCs 
and non-CCCs. 

An average temperature profile region is derived by averag- 
ing the 1-cr boundary of the scaled radial temperature profiles. 
The average temperature profile for the whole sample gives 
T(r) oc ^o.2i±o.04^ ^jjj fQj. jjjg subsample T(r) oc rO-38±(' ()4^ 

respectively, in the r < Q.lrsoo region. Note the segrega- 
tion could be more pronounced considering the spatial reso- 
lution and PFS effects. This temperature behavior in the cluster 
cores for the CCC subsample is very similar to the behavior 
for the nearby CCC sample in Sanderson et al. (2006) giving 
T(r) oc based on Chandra observations. In the outskirts 
(0.2r5oo < r < rsoo), the whole sample shows a self-similar 
behavior giving T{r) oc ^-0-36±o.i8 ^j^j^ sc^f^gj- within 20%. A 
temperature profile decreasing down to 80% of the maximum 
value with an intrinsic scatter of ~ 20% has been observed at 
about rsoo for the average of the sample. 

Studies of the cluster temperature distributions (e.g. 
Markevitch et al. 1998; De Grandi & Molendi 2002; Vikhlinin 
et al. 2005; Zhang et al. 2004a, 2006; Pratt et al. 2007) indi- 
cate a universal temperature profile with a significant decline 
beyond an isothermal center The average temperature profile 
of the pilot LoCuSS sample is consistent with the average pro- 
files from ASCA in Markevitch et al. (1998), BeppoSAX in De 
Grandi & Molendi (2002) and XMM-Newton in Zhang et al. 
(2004a, 2006) and Pratt et al. (2007), but is slighdy less steep 
than the profile from Chandra (using an assumed uncertainty of 
20% of the averaged temperature profile as an approximate il- 
lustration) in Vikhlinin et al. (2005). A similarly universal tem- 
perature profile is indicated in outskirts of clusters by simula- 
tions (e.g. Borgani et al. 2004; Borgani 2004). 



4.2. Scaled metallicity profiles 

The metallicity profiles are shown in Fig. [3] with their radii 
scaled to rsoo. We observed metallicity peaks towards the clus- 
ter centers within 0.2 rsoo for the CCCs. There is no evident 
evolutionary effect comparing the pilot LoCuSS sample to the 
nearby CCCs in De Grandi & Molendi (2002) within the mea- 
sure uncertainties. 

4.3. Scaled surface brightness profiles 

In the standard self-similar model, gas mass scales with mass 
and thus temperature as Mgas oc M oc r'-^, which gives 
5x oc r*^^ (e.g. Arnaud et al. 2002). The Sx-T scaling based 
on the empirical Mgas-T" relation provides the least scatter in 
the scaled surface brightness profiles (e.g. Neumann & Arnaud 
2001; Arnaud et al. 2002). We thus applied the current em- 
pirical relation Mg^s oc r'^ (e.g. Mohr et al. 1999; Vikhlinin 
et al. 1999; Castillo-Morales & Schindler 2003) corresponding 
to the scaling Sx T^-^ to scale the surface brightness pro- 
files (Fig.|4|i. We found a less scattered self-similar behavior at 
R > 0.2r5oo for the scaled surface brightness profiles compared 
to the profiles scaled by 5x T"-^. 

The bolometric X-ray luminosity (here we use the 0.01- 
100 keV band) is given by Lboi oc J A{r)nl{r)dV, practically 
an integral of the X-ray surface brightness to 2.5r5oo (L'^^i in 
Tabl^. The value 2.5r5oo is used because that the extrapolated 
virial radii are about 2.2-2.6r5oo. We note the luminosity only 
varies within 3% setting the truncation radius to the value from 
rsoo to 2.5r5oo. 

The surface brightness profiles show small core radii for 
the CCCs and large core radii for the non-CCCs. Similar to 
the REFLEX-DXL sample (Zhang et al. 2006) the core radii 
populate a broad range of values up to 0.2 rgoo. The X-ray lu- 
minosity is sensitive to the presence of the cool core. It can thus 
be used to probe the evolution of the cool cores. We present the 
bolometric luminosity including and excluding the R < 0.2r5oo 
region (L'^^^ and L^^") in Table[T] The fractions of the X-ray lu- 
minosity attributed to the R < Q.lrsoo region span a broad range 
up to 70%. This introduces large uncertainties in the use of the 
luminosity U^^^ as a mass indicator. The normalization varies 
with the fraction of the CCCs in the sample and the significance 
of their cool cores by a factor of 30-70%. The X-ray luminosity 
excluding the < 0.2r5oo region (^^0^^ much less dependent 
on the presence of the cool core. However, to better re-produce 
the cluster luminosity and the normalization of the luminosity- 
temperature/mass relation, we use the X-ray luminosity cor- 
rected for the cool core by assuming 5x(^) = 5x(0.2r5oo) for 
the R < 0.2r5oo region (L™;"^ in Table [T]l. Using L™J^ for the 
X-ray scaling relations we obtained reduced scatter and con- 
sistent normalization for the CCCs and non-CCCs as shown 
later. Note the luminosity is lower (by up to 10%) assuming a 
constant luminosity in the cluster core instead of a beta model, 
especially for the pronounced CCCs. However, the luminosity 
corrected for the cool core by the beta model still introduces 



^ Note this assumption was used only in the calculation for Lj 
not for any other quantity. 
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relatively significant scater dominated by the CCCs due to the 
correlation between the core radius and the slope parameter. 

4.4. Scaled cooling time profiles 

The cooling time is derived by the total energy of the gas di- 
vided by the energy loss rate 

l.SriaT 

tc = (5) 
niK 

where A, Wg, «e and T are the radiative cooling function, gas 
number density, electron number density and temperature, re- 
spectively. We compute the upper limit of the age of the clus- 
ter as an integral from the cluster redshift z up to z = 100. 
Cooling regions are those showing cooling time less than the 
upper limit of the cluster age. The boundary radius of such a 
region is called the cooling radius. The cooling radius is zero 
when the cooling time in the cluster center is larger than the 
upper limit of the cluster age. The cooling time at the resolved 
inner most radii of the surface brightness profiles and cooling 
radii are given in Table |2] 

In Fig.|5] we show the cooling time profiles with their radii 
scaled to rgoo. The CCCs show larger cooling radii in unit 
of rsoo. In the cluster centers, the cooling time is all shorter 
than the upper limit of the cluster age. The CCCs show much 
shorter cooling time than for the non-CCCs. The scaled cool- 
ing time profiles show a self-similar behavior above 0.2r5oo for 
the CCC and non-CCC subsamples, respectively. The best fit 
power law above 0.2r5()o gives tdr) oc ^i-6i±o oi for the whole 
sample. For the 4 CCCs, the best fit power law above 0.2r5oo 
gives fc(r) oc ^i-7o±o.oi p^j. j-j^g non-CCCs, the best fit power 
law above 0.2r5oo is more flat, t^^ir) oc ^i-54±o.oi rpj^g central 
cooling time is similar to the work in Bauer et al. (2005) us- 
ing the Chandra data for the same clusters towards the most 
inner bin resolved. For the pilot LoCuSS sample, the slope of 
the cooling time profile within 0.2r5oo becomes steeper with in- 
creasing significance of the cool core. This was also observed 
for the nearby clusters based on Chandra data in Sanderson et 
al. (2006) and Dunn & Fabian (2006). 

The cooling time is calculated towards the cluster center 
to the inner most bin that the surface brightness profile is re- 
solved. There can be large uncertainties in the cooling time at 
the inner most bins where the temperature measurements are 
not resolved. We note that the spatial resolution is important to 
calculate the proper cooling time in the cluster center When the 
cluster core is not well resolved, the value calculated from the 
measured temperature and electron number density gives the 
upper limit of the cooling time, specially for the CCCs show- 
ing cusped surface brightness profiles. 

4.5. Scaled entropy profiles 

The entropy is the key to the understanding of the thermody- 
namical history since it results from shock heating of the gas 
during cluster formation. The observed entropy is generally de- 
fined as S - Tn^ for clusters (e.g. Ponman et al. 1999). 
Radiative cooling can raise the entropy of the ICM (e.g Pearce 



et al. 2000) or produce a deficit below the scaling law (e.g. 
Lloyd-Davies et al. 2000). For all 12 galaxy clusters, the tem- 
perature data points are resolved below O.lrsoo. The entropy at 
0.1r2oo, -So ii^oo, can thus be used as an indicator of the central 
entropy. 

According to the standard self-similar model the entropy 
scales as 5 oc r (e.g. Ponman et al. 1999). We investigated the 
entropy-temperature relation (S-T) using 5o.ir2oo the central 
entropy and ro.a-o.srjoo the cluster temperature, ^o.irjoo indi- 
cates the physics on core scales and ro.2-o.5)5„o on global scales. 
Therefore the CCCs show significantly lower central entropies 
compared to the S-T scaling law. This can be due to that lower 
temperature systems show more pronounced cool cores cor- 
responding to lower central entropies. The entropy at O.lraoo 
can thus be used as a mechanical educt of the non-gravitational 
process which introduces large scatter of the 5o.inoo~^ relation 
(left panel of Fig. |6]l. 

Ponman et al. (2003) suggested to scale the entropy as 
S oc r"-^^ based on the observations for nearby clusters, with 
which the pilot LoCuSS sample agrees (see the left panel in 
Fig.|6]l. 0.1r2()o is slightly smaller than 0.2r5oo for the LoCuSS 
sample. To avoid the extrapolation to calculate r2oo, we also 
used the entropies at 0.2r5oo, So,2r^g„- The scatter is reduced by 
almost a factor of 2 for the S-T relation using So.iisoo instead 
of 5o.ii2oo because that the entropy profiles are more self-similar 
beyond 0.2r5oo than in the cluster cores. We show the S-T re- 
lation using S'o.irsoo in the right panel of Fig. |6] and Table |5] 
Within the error, the So,2tsoo~T relations for the pilot LoCuSS 
sample and the REFLEX-DXL sample cannot rule out the stan- 
dard self-similar model S cc T while also being consistent with 
Ponman et al. (2003). The CCCs (5o.2r5„o ^ T°7^^°-^^) agree 
better with the empirically determined scaling (Ponman et al. 
2003) S oc r" ''^ As shown in Fig. [T] the further away from 
the cluster center, the less scatter there is in the S-T relation 
(also see Pratt et al. 2006). There is no noticable evolution in 
the S-T relation (e.g. using 5o.3i5o„ and 5o.3r2oo) comparing the 
pilot LoCuSS sample to the REFLEX-DXL sample (z ~ 0.3, 
Zhang et al. 2006) and the nearby relaxed cluster sample in 
Pratt et al. (2006). With the redshift correction, the combined 
fit of the pilot LoCuSS sample, the REFLEX-DXL sample and 
the nearby relaxed cluster sample in Pratt et al. (2006) gives 
5o.3i-5o„ <^ and 5o.3r,o„ °^ r°-^4±o.os^ respectively. Note 

there is the segregation between the CCCs and non-CCCs in 
which the CCCs show low normalization of the S-T relation. 
With the combined data, at the high mass end the clusters are 
un-biased to CCCs and at the low mass end the clusters are bi- 
ased to CCCs (relaxed clusters in Pratt et al. 2006). The slope 
of the S-T relation for the combined data is thus steeper than 
the slope for the CCCs or non-CCCs. This is less significant for 
the S-T relation using the entropies at larger radii, e.g. 5o.3i,oo- 

We scaled the radial entropy profiles using the empirically 
scaling (Ponman et al. 2003) S oc r"^^ and rsoo. As shown in 
Fig.|8] the scaled entropy profiles for the pilot LoCuSS sample 
are self-similar above 0.2r5()o and show the least scatter around 
0.2-0.3r5oo. However, we found the redshift correction is re- 
quired to obtain an agreement on global radial scales for the 
pilot LoCuSS sample, Birmingham-CfA sample (Ponman et al. 
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2003), REFLEX-DXL sample (Zhang et al. 2006) and the clus- 
ter sample in Pratt et al. (2006). 

After the redshift correction the combined entropy profiles 
of the pilot LoCuSS sample give the best fit S(r) oc 
above 0.2r50(). A similar power law was found as 5 oc r*^-^^ by 
Ettori et al. (2002), S oc r" ''^ by Piffaretti et al. (2005), and 
Sir) oc r^oo±o.07 f^j. jj^g refLEX-DXL sample (> 0.2r5oo)- 
The spherical accretion shock model predicts S oc r^^ (e.g. 
Tozzi & Norman 2001; Kay 2004). The combined fit of the 
entropy profiles for the 4 CCCs in the pilot LoCuSS sample 
gives S oc ^i io±o o5^ similar to the prediction of the spherical 
accretion shock model and consistent with the relaxed nearby 
clusters in Pratt et al. (2006) giving S oc r"'^. Since the clus- 
ters appear more "relaxed" in Pratt et al (2006) the entropy 
profiles of their sample are very consistent with the CCCs in 
the pilot LoCuSS sample. The combined fit of the entropy pro- 
files for the 8 non-CCCs in the pilot LoCuSS sample gives 

S oc /WViO.OS^ 

The different slopes of the entropy profiles for the CCCs 
and non-CCCs may indicate the different phases or different 
origins of the clusters with respect to the cluster morpholo- 
gies. Observations suggest an evolution of the cluster mor- 
phology such that mergers happen more frequently in high 
redshift galaxy clusters (e.g. Jeltema et al. 2005; Hashimoto 
et al. 2006). The non-CCCs are dynamically young and thus 
show the evidence of the recent mergers judging from their X- 
ray flux images such as elongation and extended low surface 
brightness features (e.g. Abell 1763 in the pilot LoCuSS sam- 
ple). Their entropy profiles can be flattened by mergers and 
thus become less steep than the prediction of the spherical ac- 
cretion shock model. There is an evidence of AGN activities in 
some clusters (e.g. Abell 1763, Hardcastle & Sakelliou 2004), 
in which the central AGN could also affect the entropy pro- 
file. On the contrary, the CCCs appear "relaxed", and their en- 
tropy profiles agree with the prediction of the spherical accre- 
tion shock model. On the other hand, the different slopes of the 
entropy profiles of the CCCs and non-CCCs can also be due 
to geometry. The non-CCCs appear more elongated due to the 
recent mergers while the CCCs more symmetric as shown in 
the XMM-Newton flux images (Fig.[T7l also see Fiss. \B.6\W7\ 
in the electronic edition of the Journal for the whole sample). 
The flattening of the entropy profiles of the non-CCCs can also 
be a geometric effect due to the azimuthal average. 

4.6. Scaled total mass and gas mass profiles 

The mass profiles were scaled with respect to M500 and rsoo, 
respectively (Fig. |9j. We found the least scatter at radii above 
0.2r50(). In the inner regions (< 0.2r5oo), the mass profiles vary 
significantly with the cluster central dynamics (Fig. |9]l. As 
found in the strong lensing studies based on high quality HST 
data (Smith et al. 2005), the cluster cores can be very com- 
plicated with multi-clump DM halos. The peculiarities in the 
cluster cores for the individual clusters may explain the large 
scatter of the mass profiles on core scales where the ICM does 
not always trace DM. This scatter may also suggest the differ- 
ent phases of the clusters. 



Similar to the mass profiles, the scaled gas mass profiles 
appear self-similar at radii above 0.2r5oo but show less scatter 
(a few per cent) than for the scaled mass profiles. 

5. X-ray scaling relations 

To use the mass function of the cluster sample to constrain cos- 
mological parameters, it is important to calibrate the scaling 
relations between the X-ray luminosity, temperature and gravi- 
tational mass, the fundamental cluster properties including also 
velocity dispersion. Massive clusters, selected in a narrow red- 
shift range are important to constrain the normalization and to 
understand the scatter in the scaling relations. Both REFLEX- 
DXL and pilot LoCuSS samples are such flux-limited, mor- 
phology unbiased samples. Comparing the X-ray scaling rela- 
tions of such samples to samples in other narrow redshift bins 
can constrain the evolution of the X-ray scaling relations. 

The scaling relations are generally parameterized by a 
power law (Y = YqX'^'). Note the pilot LoCuSS clusters are 
in a narrow temperature range, which has the disadvantage to 
constrain the slope parameter To study the normalization seg- 
regation of the scaling relations for the CCCs and non-CCCs, 
we fit the normalization with fixed slope parameters as often 
used in previous studies for the pilot LoCuSS sample and the 
pilot LoCuSS non-CCCs. For each relation, we also performed 
the fitting with both normalization and slope free. The fitting 
slopes generally appear consistent with the fixed slopes within 
the errors. The scatter describes the dispersion between the 
observational data and the best fit. We list the best fit power 
law and the scatter in Table |5] For the pilot LoCuSS sample 
the cluster masses are uniformly determined from high quality 
XMM-Newton data. This guarantees the minimum systematic 
error due to the analysis method. 

5.1. Cool cores and the X-ray scaling relations 

In the cluster cores we observed broad scatter of the X-ray 
properties with respect to the cluster morphologies (e.g. CCCs 
and non-CCCs) for the scaled profiles of the X-ray quantities 
(e.g. temperature, metallicity, surface brightness, cooling time, 
entropy, total mass and gas mass). This is most probably due to 
the effects of different physical processes rather than statistical 
fluctuations in the measurements (Zhang et al. 2004b, 2005b; 
Finoguenov et al. 2005). The CCCs and non-CCCs can be the 
clusters at their different phases. 

The X-ray luminosity integrated over the whole cluster re- 
gion introduces significant normalization segregation between 
CCCs and non-CCCs, and large scatter of the X-ray scaling 
relations dominated by non-CCCs. It is thus worth to take a 
closer look how the cool cores affect the scaling relations and 
their scatter as follows. 

For SoAr^og-T relation, including the CCCs in the sample 
introduces not only a steeper slope but also a lower normaliza- 
tion (~ 30%). 

Using the X-ray luminosity including the cool core (L'"'^'^), 
the normalization of the L-T and L-M relations excluding the 
4 CCCs is reduced by 40%. If the temperature is measured 
also including the cool core, the normalization segregation will 
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be even larger (by a few per cent) for the L-T relation. Using 
the luminosity corrected for the < 0.2rgoo region and temper- 
ature excluding the < 0.2r5oo region, the normalization of the 
scaling relations agrees to better than 10% for the CCCs and 
non-CCCs. In such a way the scaling relations among X-ray 
luminosity, cluster mass and temperature, and their scatter are 
insensitive to the exclusion of CCCs for the pilot LoCuSS sam- 
ple. The above results show the following picture. CCCs usu- 
ally appear more "relaxed", and are thus considered to provide 
reliable X-ray mass measurements (e.g. Pierpaoli et al. 2001; 
Allen et al. 2004). At the same time, the properties of the cool 
cores exhibit the largest scatter. Therefore, the CCCs were usu- 
ally selected and the cool cores were often excised or corrected 
to study the X-ray scaling relations. 

With the bolometric X-ray luminosity including and ex- 
cluding the < 0.2rgo() region (Table|2]i, we show the normalized 
cumulative cluster number count as a function of the fraction 
of the luminosity attributed by the cluster core (< 0.2r5oo) in 
Fig. [To] Up to 70% of the bolometric X-ray luminosity is con- 
tributed by the cluster core (< 0.2rg(K)). The evolution of the 
cool cores is not significant as also explained using the temper- 
ature drop towards the cluster center in Sect. 13.21 

As mentioned above the global temperature was calculated 
excluding the < 0.2r5()o region and the X-ray luminosity was 
corrected for the < 0.2r50() region, respectively. In such a way 
the scatter in the X-ray scaling relations is minimized (< 20%) 
and the normalization of the X-ray scaling relations is insensi- 
tive to the exclusion of the CCCs (< 10%). 

5.2. Mass-temperature relation 

Given the virial theorem (T <x Mgoo/rioo) and the spherical 
collapse model (Mgoo cc r^Q^) one obtains M500 ''c r'-^. This 
scaling was also observed for nearby clusters (e.g. Finoguenov 
et al. 2001; Reiprich & Bohringer 2002; Arnaud et al. 2005). 
To compare the pilot LoCuSS sample to the existing nearby and 
more distant samples, we present the Msqq-T relation fixing the 
slope parameter to the often used value (1.5) in Fig. [TT] The 
scatter in the cluster mass in the M^oq-T relation is within 20% 
for the pilot LoCuSS sample. 

Evrard et al. (1996) simulated ROSAT observations of 58 
nearby clusters [z ~ 0.04, 1-10 keV), for which the normal- 
ization of the M-T relation agrees with the pilot LoCuSS sam- 
ple. Chen et al (2007) investigated ROSAT and ASCA obser- 
vations of a flux-limited, morphology-unbiased sample of 106 
nearby clusters (z < 0.15, 1-15 keV, extended HIFLUGCS). 
The scaling relation of the massive non-CCCs (T > 3 keV) 
in their sample shows an excellent agreement with the pilot 
LoCuSS sample. This also indicates that the cool core correc- 
tion for the pilot LoCuSS sample introduces less bias into the 
normalization towards the CCCs. Note the early work related 
to the HIFLUGCS sample (Finoguenov et al. 2001; Popesso et 
al. 2005) gives similar scaling relations as in Chen et al. (2007) 
for the whole HIFLUGCS sample. Vikhlinin et al (2006a) de- 
rived the M-T relation for 13 low-redshift clusters (z < 0.23, 
0.7-9 keV) using Chandra observations, which normalization 
agrees with ours within the scatter (20%) for the pilot LoCuSS 



sample. The scaling in Arnaud et al. (2005) is based on XMM- 
Newton observations of 6 relaxed nearby clusters (z < 0.15, 2- 
9 keV, T > 3.5 keV). The temperatures used in their work were 
determined in the 0.1-0.5r2oo region, which also excludes the 
cluster cores as we did using the 0.2-0.5r5oo region. The tem- 
perature profiles in Arnaud et al. (2005) are almost flat in the 
0.5r5()o-0.5r2oo regions and their values are close to the max- 
imum temperatures. The measured global temperature is thus 
higher when the 0.5r5oo-0.5r20{) region is used as in Arnaud et 
al. (2005) than the global temperature measured excluding the 
0.5r5oo-0.5r2oo- We, also in Zhang et al. (2006), measured the 
global temperatures using the region only up to 0.5rg(K) deter- 
mined by the photon statistic. Therefore the normalization of 
the M-T scaling relations in this work and Zhang et al. (2006) 
is slightly higher than in Arnaud et al. (2005). 

The M-T relation in Ettori et al. (2004) is based on 
Chandra observations of 28 high redshift clusters (0.4 < z < 
1.3, 3-11 keV). It shows a higher normalization compared to 
the above published relations we mentioned. This could be due 
to the different way to define the density contrast used in the 
redshift evolution correction and the method to calculate the 
global temperatures in Ettori et al. (2004). We used the redshift 
dependent density contrast (Bryan & Norman 1998) in the red- 
shift evolution correction. We note that in this work the cool 
cores are excluded in calculating the global temperatures. In 
Ettori et al. (2004) the central regions including cluster cool 
cores were used to measure the global temperatures, which 
tends to give lower temperatures for given cluster masses for 
the clusters showing cool cores. Therefore the normalization 
of the M-T relation in Ettori et al. (2004) can be higher. The 
scaling relation in Zhang et al. (2006) is based a flux-limited 
morphology-unbiased sample of 13 medium distant clusters 
(0.26 < z < 0.31, r > 5 keV, REFLEX-DXL) observed 
by XMM-Newton. Compared to the sample in this work the 
REFLEX-DXL sample shows a slightly higher normalization 
but well within the mass scatter (~ 20%) for the pilot LoCuSS 
sample. 

No evident evolution is found for the M-T relation com- 
paring the pilot LoCuSS sample to the nearby and more distant 
samples within the scatter. 

5.3. Gas mass-temperature relation 

Assuming the gravitational effect dominates in galaxy clus- 
ters, the gas should follow the collapse of the DM giving 
Mgas.soo M5Q0 and thus Mgas.soo 0^ r'-^ (e.g. Arnaud 2005). 
However, the non-gravitational effects become important for 
the ICM such that the shape of the ICM density profile de- 
pends on the cluster temperature (e.g. Ponman et al. 1999) and 
the observed slope of the Mgas.soo-^ relation becomes steep 
(~ 1.8, e.g. Mohr et al. 1999; Vikhlinin et al. 1999). The slope 
of the gas mass-temperature relation for the pilot LoCuSS sam- 
ple is consistent with the previously observed value 1.8 within 
the eiTor. Therefore, we present the Mgas.soo-^ relation for the 
pilot LoCuSS sample with the fixed slope of 1.8 to compare 
with the recently published results in Fig.[TT] The scatter of the 
gas mass for the Mgas.soo-^ relation is within 15% for the pilot 
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LoCuSS sample. Note the weighted (unweighted) fitting slope 
of the Mgas,5oo-r relation is 2.8 + 0.8 (1 .6 + 0.2). 

The Mgas,5oo-7^ relation in Mohr et al. (1999) is based on 
a flux-limited sample of 45 nearby clusters spanning a tem- 
perature range of 2-10 keV observed by ROSAT and ASCA. 
Castillo-Morales & Schindler (2003) investigated a sample of 
10 nearby clusters (0.03 < z < 0.09, 4.7-9.4 keV) also ob- 
served by ROSAT and ASCA. The Mg^.^^oo-T relation for the 
pilot LoCuSS sample agrees with the published relations for (i) 
the nearby samples in e.g. Mohr et al. (1999), Castillo-Morales 
& Schindler (2003), and Chen et al. (2007, the non-CCCs with 
r > 3 keV in HIFLUGCS), and (ii) the more distant sample in 
Zhang et al. (2006, REFLEX-DXL). 

5.4. Luminosity-temperature relation 

In the standard self-similar model Mgas,500 °^ A^soo oc r'^ and 
M500 oc r^QQ give L oc Mg^^r^-^/r^Qg oc T^. However, the L- 
T relation deviates from the L oc scaling due to the de- 
pendence of the ICM distribution on the cluster temperature 
(e.g. Neumann & Arnaud 2001) and becomes steeper. Given 
the empirical scaling Mgas.soo °^ T^-^, one gets L oc T^-^ as of- 
ten found in observations. We thus fixed the slope to the often 
observed values (e.g. 2.6 for ioT2 4keV-^ 2.98 for L^- 
T, e.g. observed in Reiprich & Bohringer 2002), and compared 
the findings for the pilot LoCuSS sample to the recent results 
in Fig. [12] For the pilot LoCuSS sample, the scatter of the lu- 
minosity for the L-T relations is within 15%. 

The L-T relation has been intensively studied for the 
nearby cluster samples using the luminosity based on ROSAT 
observations and the temperature based on ASCA observa- 
tions, for example, in (i) Markevitch (1998, 30 clusters, 0.04 < 
z < 0.09, r > 3.5 keV), (ii) Evrard & Arnaud (1999, 24 
clusters, z < 0.37, T > 2 keV), (iii) Reiprich & Bohringer 
(2002, in which the flux-limited morphology-unbiased sample 
HIFLUGCS was initially constructed and investigated, note the 
fits are for the whole HIFLUGCS sample including groups), 
(iv) Ikebe et al. (2002, a flux-limited sample of 62 clusters, 
z < 0.16, 1-10 keV), and (v) Chen et al. (2007, HIFLUGCS). 

The L-T relation deviates from the standard self-similar 
prediction L oc T^, but shows no evolution comparing the 
pilot LoCuSS sample to the nearby samples in Markevitch 
(1998) and Arnaud & Evrard (1999) using representative 
non-CCCs. As flux-limited morphology-unbiased samples, the 
Lg°[^2 4kev~^ relation for the pilot LoCuSS sample agrees with 
the HIFLUGCS sample (e.g. Reiprich & Bohringer 2002, in- 
cluding groups; Chen et al. 2007, non-CCCs with T > 3 keV), 
the sample in Ikebe et al. (2002) and the REFLEX-DXL sam- 
ple (e.g. Zhang et al. 2006). Therefore no evident evolution is 
observed for the L-T relation up to redshift 0.3. 

Kotov & Vikhlinin (2005) applied an alternative redshift 
evolution for the L-T relation of 10 distant non-CCCs (0.4 < 
z < 0.7, T > 3.5 keV), giving L oc t^M(i + j;)i-8. We applied 
the redshift coiTection to the 7 clusters available in Kotov & 
Vikhlinin (2005) also observed by XMM-Newton and found 
an agreement with the pilot LoCuSS sample within the scat- 
ter. We note that in Kotov & Vikhhnin (2005) the tempera- 



ture was measured in the 70 kpc-rsoo region and luminosity 
in the 70-1400 kpc region, in which a small central region was 
excluded in calculating the temperature/luminosity. This tends 
to give slightly lower/higher temperature/luminosity and thus 
higher normalization of the L-T relation compared to the pilot 
LoCuSS sample. On the other hand, Kotov & Vikhlinin (2005) 
used the same aperture for the correction which gives a more 
significant correction of the cool core for a less massive system 
and could thus provide lower normalization of the L-T relation 
for less massive systems. 

5.5. Luminosity-mass relation 

Witii M oc T^-^ and L oc T^ in the standard self-similar model 
one derives L oc M^-^^. With the observed scaling relations 
in e.g. Reiprich & Bohringer (2002), Lo.i-2.4keV °^ T^'^ gives 
Lo.i-2.4keV M'-^^ and Lboi T^-^^ gives Lboi M^-^^, re- 
spectively. The fitting of the pilot sample favors the observed 
scaling relations. Therefore, we fixed the slope to 1.73 for 
^oT-2.4keV-^ and 1.99 for L^^-T, and compared the pilot 
LoCuSS sample to the recent results in Fig.[T3] The scatter of 
the luminosity and mass for the L-T relations is within 15% 
and 20%, respectively. 

The pilot LoCuSS sample agrees with the nearby sample 
HIFLUGCS (Reiprich & Bohringer 2002; Popesso et al. 2005; 
Chen et al. 2007, non-CCCs with T > 3 keV) and the more 
distant sample REFLEX-DXL (Zhang et al. 2006) also as flux- 
limited morphology-unbiased samples. Therefore we observed 
no evident evolution of the luminosity-mass relation. 

5.6. rjoo vs. global temperature 

As found in Evrard et al. (1996) the standard self-similar 
model predicts rsoo 0^ M'^^ oc r"^. The best fit for the pilot 
LoCuSS sample confirms the standard self-similar model, and 

r,U,oo ^ - 1 A-0.335±0.1{)6x0.57±().12 e'/,\-I/-a /a \-0.5 A/r„„ 
gives rsoo - 10 0.2-().5;-5oo (^c,z/Ac.,o) MpC. 

It agrees with the scaling relation, rsoo oc ^O-SOiO Os ^^j. g j.g_ 
laxed nearby clusters (z < 0.15, 2-9 keV, T > 3.5 keV) in 
Arnaud et al. (2005) also observed by XMM-Newton. 

The above comparison of the X-ray scaling relations in this 
work with the published results (Figs. fTTHlJl l for the nearby 
and more distant samples shows that the evolution of the scal- 
ing relations can be accounted for by the redshift evolution 
given in Sect. |4] As shown in Table|5]the normalization of the 
scaling relations agrees for the CCCs and non-CCCs to better 
than 10% using the temperature measured excluding the cool 
core and the luminosity coiTected for the cool core. The X-ray 
quantities such as temperature and luminosity can be used as 
reliable indicators of the cluster mass with the scatter less than 
20%. In general, the slopes of the scaling relations indicate the 
need for non-gravitational processes. This fits into the general 
opinion that galaxy clusters show a modified self-similarity up 
to z ~ 1 (e.g. Arnaud 2005). As also demonstrated in simula- 
tions (e.g. Poole et al. 2007), mergers only alter the structure 
of compact cool cores, and the outer structure (> 0.1r2oo or 
> 0.2r5oo) is survived after the merger events. The scaling rela- 
tions are thus preserved for galaxy clusters. 
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6. Discussion 

6.1. Comparison of the X-ray results 

The Chandra data of Abell 1689 (Xue & Wu 2002), Abell 383 
and Abell 2390 (Vikhlinin et al. 2006a) were analyzed in de- 
tails. Bauer et al. (2005) present the cooling time of 8 clusters 
in the pilot LoCuSS sample based on Chandra data. Smith et al. 
(2005) measured the Chandra temperatures for 9 clusters in this 
sample. The published Chandra results for these clusters agree 
with the XMM-Newton results in this work. We also obtained 
an agreement with the published results based on the same 
XMM-Newton data, e.g. Abell 1689 (Andersson & Madejski 
2004), Abell 1835 (Majerowicz et al. 2002; Jia et al. 2004) and 
Abell 2218 (Pratt et al. 2005). 

6.2. Gas profiles in the outskirts 

The generally adopted y6 model (f5 = 2/3) gives rig oc r"^. 
However, Vikhlinin et al. (1999) found a mild trend for p 
to increase as a function of cluster temperature, which gives 
y6 ~ 0.80 and /ig °^ r^^''^ for clusters around 10 keV. Bahcall 
(1999) also found that the electron number density scales as 
rif. oc r"^ "* at large radii. Zhang et al. (2006) confirmed their 
conclusion that «e °^ r^^-'^^ for the REFLEX-DXL clusters. Due 
to the gradual change in the slope, one should be cautious to use 
a single slope double-/? model which might introduce a system- 
atic eiTor in the cluster mass measurements (as also described 
in e.g. Horner 2001). Similarly, we performed the power law fit 
of the ICM density distributions at radii above 3' for the present 
sample and obtained an average of oc ^-2.2±o.i 

6.3. Lensing to X-ray mass ratios 

Smith et al. (2005) present a uniform strong lensing analysis of 
9 clusters in this sample. We took these 9 overlapping clusters 
as a subsample (hereafter the S05 subsample. X-ray selected) 
and compared the strong lensing and X-ray masses. Ten clus- 
ters in the pilot LoCuSS sample have CFH12k data. The de- 
tailed data reduction method can be found in Bardeau et al. 
(2005) and the weak lensing results of the sample in Bardeau 
et al. (2006), respectively. The CFH12k optical luminosity 
weighted galaxy number density contours and weak lensing 
masses are given in Bardeau et al. (2006). We took these 10 
clusters as a subsample (hereafter the B06 subsample. X-ray 
selected) and compared the weak lensing and X-ray masses. In 
total eleven clusters have both X-ray and weak/strong lensing 
masses. 

To compare to the lensing mass we calculated the projected 
X-ray mass distribution for each cluster using its observed ra- 
dial mass profile. We note that the clusters are detected up to the 
radii between rsoo and r20(). Because the extrapolated virial radii 
are about 2.2-2.6r50(), we projected the observational mass pro- 
file using a truncation radius of 2.5r5oo. We also performed the 
projection to the weak lensing determined r^Qp. The comparison 
of the projected X-ray masses using these 2 different truncation 
radii, 2.5r5oo and ?"Jqq), shows that the X-ray projected mass is 
insensitive to the truncation radius, varying within a few per 



cent. The eiTor introduced by the projection to 2.5r5oo is rela- 
tively small compared to the projected mass uncertainties, for 
example, ~ 30% at rJo^, for Abell 68 using the uncertainties 
of the electron number density and temperature measurements 
by Monte Carlo simulations. Therefore we adopt the projected 
X-ray mass using the truncation radius of 2.5r50() for further ap- 
plications. We compare the X-ray masses to the strong lensing 
masses in Smith et al. (2005) at the X-ray determined rasoo and 
to the weak lensing masses in Bardeau et al. (2006) at the weak 
lensing determined r^^, respectively (Fig. [14] and Fig.fTsll. The 
scatter in the strong/weak lensing to X-ray mass ratio is com- 
parable. 

Within the error, 6 out of 9 clusters show consistent mass 
estimates between X-ray and strong lensing approaches. One 
third of the clusters show higher strong lensing masses com- 
pared to the X-ray masses up to a factor of 2.5. The mean of 
the strong lensing to X-ray mass ratio is 1.53 with its scatter 
of 1.07. On average the strong lensing approach tends to give 
larger masses than the X-ray approach especially for dynami- 
cally active clusters. 

Within the error, 4 out of 10 clusters show consistent mass 
estimates between X-ray and weak lensing approaches. The 
mean of the weak lensing to X-ray mass ratio is 1.12 with 
its scatter of 0.66. Among the clusters showing inconsistent 
masses, half have higher weak lensing masses compared to 
the X-ray masses. On average the X-ray and weak lensing ap- 
proaches agree better. 

6.4. Mass discrepancy 

The mass discrepancy between X-ray and lensing can be a 
combination of the measurement uncertainties and the physics 
in the individual clusters. The following issues have to be con- 
sidered in properly comparing the lensing and X-ray masses. 
(1) Simulations in Poole et al (2006) show that the tempera- 
ture fluctuations AT/T ~ 20% can pass the virialization point 
and persist in the compact cluster cores. Therefore tempera- 
ture fluctuations at the 20% level do not necessarily indicate 
a disturbed system, but introduce uncertainties in the cluster 
mass estimates in the central region. (2) As demonstrated in 
Puy et al. (2003), the relative error of the X-ray surface bright- 
ness goes almost linearly as a function of the rotating angle 
and axis ratio of the cluster, respectively. The projection of 
the X-ray masses here is based on the assumption of spheri- 
cal symmetry, in which the total mass can be easily overesti- 
mated/underestimated for aspherical clusters. (3) Specially for 
the mass distribution in a CCC the cool core has to be well 
resolved to guarantee the proper X-ray mass estimate in the 
central region. (4) The line-of-sight merger can enhance the 
possibility to observe the cluster by strong lensing. A strong 
lensing selected sample can be boosted because of this. (5) 
The lack of reUable redshift measurements of the faint galax- 
ies in the weak lensing analysis causes confusion of the lensed 
background galaxies and unlensed cluster galaxies. This intro- 
duces a significant uncertainty in the weak lensing mass esti- 
mate (Pedersen et al. 2006). (6) The lack of reliable arc redshift 
measurements causes uncertainties in the strong lensing mass 



Zhang et al.: Scaling Relations and Mass Calibration of the X-ray Luminous Galaxy Clusters at z=0.2 



13 



estimate. (7) The mass distribution of the DM halo can be more 
complex with multi-components which can significantly affect 
the strong lensing analysis (e.g. Smith et al. 2005). 

The mass discrepancy between X-ray and lensing partially 
come from the interpretation of galaxy dynamics and ICM 
structure in galaxy clusters (e.g. Smith et al. 2005; Zhang et al. 
2005a). The optical cluster morphology is described in Bardeau 
et al. (2006). To reveal the link between the X-ray lensing mass 
discrepancy and the characters of the individual clusters, we 
investigated the strong/weak lensing to X-ray mass ratio as a 
function of the X-ray properties such as temperature, metallic- 
ity, central entropy, cooling radius, cooling time and X-ray lu- 
minosity. However, no evident correlations are observed. This 
could be a result of the low statistic of only 10 clusters in the 
subsample. A large sample with better statistic could be helpful 
for better understanding. In Fig.[T6]we present the projected gas 
mass fractions at r25oo using the strong lensing masses and at 
rjQg using the weak lensing masses, which span a broad range 
up to 0.3. The scatter of the projected gas mass fractions using 
the strong (weak) lensing approach is larger by a factor of 1 .5 
(2.5) than the scatter using the X-ray masses. 

Combining data on the optical cluster morphology and the 
X-ray appearance of the cluster may tell more about the cluster 
structure and mass discrepancy. The soft band images are less 
temperature dependent and thus better reflect the ICM density 
distribution. The MOSl data show less serious gaps which is 
important to check the cluster morphology. To have a closer 
look at the cluster morphology we created flat fielded point 
source included XMM-Newton MOS 1 flux images in the 0.7- 
2 keV band binned in 8" x8". The weighted Voronoi tesselation 
method (Cappellari & Copin 2003; Diehl & Statler 2006) was 
used to bin the image to 5cr significance. In Fig. [T7](aZio see 
Figs. \B.6HB.7\ in the electronic edition of the Journal for the 
whole sample), we superposed the CFT[12k optical luminosity 
weighted galaxy density contours (Bardeau et al. 2006) on the 
X-ray adaptively binned images. 

6.4.1 . Strong lensing and X-ray masses 

The morphology of the X-ray images overlaid with optical con- 
tours can be used to understand the mass discrepancy in the 
central regions combined with the substructure fraction in the 
strong lensing analysis (Smith et al. 2005). We thus present the 
strong/weak lensing to X-ray mass ratio versus the substruc- 
ture fraction (Fig.fTSll. Three clusters show discrepancies of the 
masses using the strong lensing and X-ray approaches. Here we 
discuss these 3 clusters, which show a significant link between 
the substructure fraction and strong lensing to X-ray mass ratio. 

Abell 383 has very low substructure fraction, appearing 
"relaxed" in the strong lensing analysis. The X-ray tempera- 
ture profile show a strong decrease towards the cluster center, 
which can be even more significant considering the PSF cor- 
rection. This can add uncertainties in the X-ray mass estimate. 

Abell 773 and Abell 2218 have very high substructure frac- 
tions showing multi-modal DM halos as found in Smith et al. 
(2005). In X-rays the multi-clumps are smeared out and both 
clusters appear symmetric. This may explain the large mass 



discrepancy between strong lensing and X-ray for Abell 773 
(by a factor of 2.3) and Abefl 2218 (by a factor of 3). Refiable 
arc redshift measurements have been obtained in the strong 
lensing mass estimate for Abell 2218 (e.g. Smith et al. 2005 
and the references in). The mass distribution of the DM halo in 
Abell 2218 shows a complex structure with multi-components 
in the central region (Smith et al. 2005). The pronounced mass 
discrepancy between X-ray and strong lensing for Abell 2218 
is mostly due to the combination of (1), (2), (4) and (7) as partly 
indicated in Girardi et al. (1997) and Pratt et al. (2005). 

6.4.2. Weak lensing and X-ray masses 

Among the 6 clusters with inconsistent weak lensing and X-ray 
masses, the mass discrepancy is not significant for Abell 68 and 
Abell 1763. 

The highest projected gas mass fractions were observed in 
Abell 963 (~ 0.3) and Abell 267 (~ 0.2) using the weak lens- 
ing masses (Fig. [T6l ). The X-ray masses are likely reliable for 
these 2 clusters since both clusters appear regular in their X-ray 
flux images. Therefore the high projected gas mass fractions 
using the weak lensing masses may indicate that the weak lens- 
ing masses are underestimated for these 2 clusters. This could 
explain the low weak lensing to X-ray mass ratio (~ 0.4) for 
AbeU 267 and Abell 963. 

Abell 1835 shows the most significant mass discrepancy 
(by a factor of 2.5) between X-ray and weak lensing. At 
'"200"^ = 2.1 1 Mpc given by the weak lensing analysis in Clowe 
& Schneider (2002), our X-ray mass is (10.4 + 2.1) x W^Mq, 
in a good agreement with their weak lensing mass (8.5 + 0.8) x 
10'"* Mq. The mass profile within 0.5r5()o agrees with the strong 
lensing mass profile in Smith et al. (2005). The X-ray mass 
(A^soo) of Abell 1835 also agrees with the previous X-ray re- 
sults in Majerowicz et al. (2002). The pronounced mass dis- 
crepancy between X-ray and weak lensing using CFH12k data 
for Abell 1835 is mostly due to the combination of (1), (2), (3) 
and (5). The mass discrepancy for Abell 383 appears similar to 
Abell 1835. 

As shown in Fig. [15] there is a trend that the weak lens- 
ing to X-ray mass ratio increases with the characteristic cluster 
size (e.g. /"joq)- The disagreement between the weak lensing 
and X-ray masses becomes worse with increasing radius. This 
could be understood as follows, (i) There is little X-ray emis- 
sion in the outskirts beyond rsoo. The extrapolation of the X-ray 
mass profiles are too steep compared to the cluster potential in 
the outskirts. The projected X-ray mass at rjgg could thus be 
underestimated, (ii) The LSS component becomes important 
at the boundary of the cluster (e.g. the REFLEX-DXL cluster 
RXCJ0014.3-3022, Braglia et al. 2007) which can enhance the 
weak lensing mass by including filaments in the projected clus- 
ter mass (Pedersen & Dahle 2006). (iii) Two colors were used 
to distinguish the lensed background galaxies and un-lensed 
cluster galaxies. The background confusion could also intro- 
duce some uncertainties, which requires reliable photometric 
redshift measurements to improve the situation (e.g. Gavazzi 
& Soucail 2007). 
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6.5. Luminosity-mass relation using lensing masses 

The purpose to calibrate the L-M relation is to seek the repre- 
sentative L-M scaling and to understand the scatter so that the 
global luminosity can be used as a cluster mass indicator. The 
global luminosity is thus used in the L-M relations, but with 
different independent cluster mass estimates. In such a way, 
the difference of the scater in the L-M relations using different 
mass approaches can help us better understanding the scatter 
and thus sources of systematics. The calibration of the L-M 
relation could thus provide the proper scatter and sources of 
systematic errors which should be included in the cluster lumi- 
nosity function for the cluster cosmology. 

We investigated the calibration of the luminosity-mass 
scaling relation using independent approaches such as X-ray 
and gravitational lensing for an X-ray selected flux-limited 
sample in a narrow redshift bin (z ~ 0.2). The X-ray clus- 
ter mass, based on the X-ray selection criteria and calculated 
from the X-ray quantities, is correlated with the X-ray luminos- 
ity. The gravitational lensing approach provides a unique tool 
to measure the cluster mass without assuming hydro-statistic 
equilibrium in the ICM required in the X-ray approach. The 
lensing mass is independent of the X-ray properties such as 
temperature, luminosity etc. Using the lensing masses instead 
of the X-ray masses in the scaling relations can guarantee less 
systematics and eiTors. The lensing and X-ray approaches can 
thus be combined to understand the scaling relations and to re- 
veal the physics in the scatter of the scaling relations. 

Smith et al. (2005) performed such studies of the mass- 
temperature relation on r250{) scales combining the Chandra and 
HST data. Petersen et al. (2006) performed similar studies for a 
collected sample from archive. With the pilot LoCuSS sample, 
there are 3 advantages to the previous studies. (1) Lensing sig- 
nals are sensitive to the two-dimensional mass distribution and 
thus tend to pick up the random structures in projection along 
the line-of-sight. The first advantage is that the subsample of 
the overlapping clusters for combining the X-ray and lensing 
studies is X-ray selected, no bias towards the clusters show- 
ing line-of-sight mergers. (2) Both the strong/weak lensing and 
X-ray data were uniformly observed and consistently analyzed 
for the overlapping clusters. The second advantage is that the 
systematics of the sample are better controlled than the sample 
collected in the archive. (3) Both strong lensing and weak lens- 
ing masses are combined with the high quality XMM-Newton 
data covering the scales extending to radii larger than rgoo. 

In Fig. [18] we show the luminosity-mass scaling relations 
using X-ray and lensing masses. We observed the correlation 
between the strong/weak lensing mass and X-ray luminosity. 
The scatter using the lensing masses is about 40%, larger than 
using the X-ray masses (< 20%). The large scatter using the 
strong/weak lensing masses can be due to the combination of 
(1) the dependence between the X-ray luminosity and X-ray 
mass since the X-ray luminosity is measured within the scale 
2. 5 rgoo determined by the X-ray approach instead of the scale 
Tjoo determined by the lensing approach, (2) the small size of 
the sample, (3) the possible additional uncertainties in the X- 
ray and lensing mass estimates due to the multi-component DM 
halo, (4) the use of a definition of the global X-ray luminosity 



such as to reduce the scatter in the X-ray scaling relations, and 
perhaps (5) the physics for the individual clusters such as cen- 
tral AGNs. 

We went through the radio archive, and found that there is 
potentially one central AGN in each cluster for the LoCuSS 
sample. There is an evidence of the link between their AGN 
activities and the scatter of the scaling relation. Abell 1763 
shows one twin-jet radio source in the center (e.g. Hardcastle 
& Sakelliou 2004), and one twin-lobe radio source ~ 5' off- 
center together with an extended low surface brightness X-ray 
feature. The archive SDSsQ spectrum shows that the cD galaxy 
of Abell 1835 is a LINER. As 2 clusters with the most signif- 
icant indication of AGN activities in the sample, Abell 1763 
and Abell 1835 coincidentally introduce the most pronounced 
scatter in the L-M scaling relations using both strong lensing 
and weak lensing masses. 

7. Summary and conclusions 

We performed a systematic analysis to measure the X-ray 
quantities based on XMM-Newton observations for the pilot 
LoCuSS sample, a flux-limited morphology-unbiased sample 
at z ~ 0.2 consisting of 12 X-ray luminous galaxy clusters. We 
investigated various X-ray properties, the scaling relations and 
their scatter, and compared the X-ray and gravitational lensing 
masses for this sample. We summarize the main conclusions as 
follows. 

(i) Self-similarity of the scaled profiles of the X-ray prop- 
erties 

An almost self-similar behavior of the scaled profiles of X- 
ray properties, such as temperature, surface brightness, cooling 
time, entropy, gravitational mass and gas mass, has been found 
for radii above 0.2 rgoo for the sample. 

- Based on XMM-Newton observations, we obtained an av- 
erage temperature profile of the sample with a plateau at 
0.2-0.5 rgoo and a drop to 80% of the maximum temper- 
ature at 0.5rgoo with 20% scatter. For the 4 CCCs, we 
observed cool gas showing lower temperatures than 70% 
of the maximum temperatures in the cluster cores. In the 
r < 0.2rgoo region, the average temperature profile for the 
whole sample gives T{r) oc ^0 2i±o.04^ ^^j. CCC 
subsample T{r) oc ^0-38±o.04^ respectively. In the outskirts 
(0.2rgoo < r < rgoo), it gives T(r) oc ^-('■36±o.i8 f^j. ^j^g ^jjole 
sample, without an evident difference between the CCCs 
and non-CCCs at Icr significance level. The averaged tem- 
perature profile for the pilot LoCuSS sample at z ~ 0.2 is 
consistent with the nearby and more distant cluster samples 
within the observational dispersion 

- We determined the XMM-Newton surface brightness pro- 
files up to at least rgoo. The surface brightness profiles of 
the non-CCCs show flat cores populating a broad range 
of values up to 0.2 rgoo. For the CCCs, the cluster cores 
are resolved with current XMM-Newton data. The surface 

http://cas.sdss.org/dr4/en/tools/explore/obj . asp 
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brightness profiles are quite self-similar R > 0.2r5oo for 
the sample. 

- The scaled cooling time profiles show an almost self- 
similar behavior above 0.2r5oo- The best fit power law 
above 0.2r5oo gives tdr) oc ^i-6i±(' ('i for the whole sample, 

tc(r) oc rl ™±0.01 foj. jjjg 4 QQQ^ fj^) ^ ^1.54±0.01 fQj. 

the 8 non-CCCs. 

- We found an empirical scaling, SQ,2rsm 7"''-8i±o.i9 jj^g ^_ 
T relation for this sample agrees with the REFLEX-DXL 
sample and the sample in Pratt et al. (2006). After the red- 
shift evolution correction, the entropy profiles for the pilot 
LoCuSS sample agree with the nearby clusters in Ponman 
et al. (2003) and Pratt et al. (2006) and the more distant 
clusters in Zhang et al. (2006) within the observational dis- 
persion. The entropy profiles at r > 0.2r5oo for the whole 
sample give S{r) oc ^i oi±o o4 n^^^ non-CCCs show a com- 
bined entropy profile of S{r) oc /) 97±o.05 rpj^g entropy pro- 
files for the 4 CCCs give, S(r) oc i()±0 05^ consistent with 
the spherical accretion shock model prediction. 

(ii) X-ray scaling relations 

- The scaling relations are sensitive to the exclusion of CCCs 
when the X-ray luminosity and temperature are measured 
including the cluster cores. The cluster cores (< 0.2r5oo) 
contribute up to 70% of the bolometric X-ray luminosity. 
Using the X-ray luminosity coiTected for the cool core and 
temperature excluding the cool core, one not only mini- 
mizes the scatter in the scaling relations but also obtains 
an agreement better than 10% between the normalization 
for the CCCs and non-CCCs. 

- For the pilot LoCuSS sample the mass scatter is less than 
20% in the mass observable scaling relations. The X-ray 
scaling relations show no evident evolution comparing the 
pilot LoCuSS sample to the nearby and more distant sam- 
ples within the observational dispersion after the redshift 
evolution correction. This fits the general opinion (e.g. 
Maughan et al. 2003; Arnaud 2005; Arnaud et al. 2005; 
Vikhlinin et al. 2006a; Zhang et al. 2006) that the evolu- 
tion of galaxy clusters up to z ~ 1 is well described by 
a self-similar model for massive clusters. With the current 
observations, the X-ray quantities such as temperature and 
luminosity can be used as reliable indicators of the cluster 
mass within 20% scatter. 

(iii) Mass calibration and the luminosity-mass relation 

The discrepancies remain between the X-ray measured 
cluster masses and the gravitational lensing masses. This can be 
due to a combination of the measurement uncertainties and the 
physics in the individual clusters. The scatter in the strong lens- 
ing to X-ray mass ratios and weak lensing to X-ray mass ratios, 
respectively, is similar. We observed the coiTelation between 
the X-ray luminosity and lensing mass. The scatter using the 
gravitational lensing masses is significant (~ 40%). This can 
be either due to the unknown physics or due to that the sample 



is too small. A large X-ray selected sample with high quality 
X-ray and gravitational lensing observations should shed light 
on it. 
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Table 1. Primary parameters. Column (1): cluster name; Col. (2): optical redshift (e.g. Smith et al. 2005; Bardeau et al. 2006); 
Cols. (3,4): sky coordinates in epoch J2000 of the cluster center; Col. (5): hydrogen column density (Dickey & Lockman 1990); 
Col. (6): truncation radius corresponding to a 5 /A^ of 3 of the observational surface brightness profile; Cols. (7-9): bolometric 
luminosity including the < 0.2r5oo region, excluding the < 0.2r5oo region, and corrected for the < 0.2r5oo region, respectively. 
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0.56 ± 0.07 


0.75 ± 0.09 


Abell 1689 


0.184 


13 


11 


29.330 


-01 


20 


26.66 


0.0182 


8.0 


2.94 ±0.10 


0.96 ± 0.06 


1.49 ±0.10 


Abell 1763 


0.228 


13 


35 


18.115 


+41 


00 


03.89 


0.0936 


7.5 


1.75 ±0.15 


1.21 ±0.12 


1.50 ±0.15 


Abell 1835 


0.253 


14 


01 


01.865 


+02 


52 


35.48 


0.0232 


7.0 


5.18 ±0.17 


1.56 ±0.10 


2.29 ±0.17 


Abell 2218 


0.176 


16 


35 


53.775 


+66 


12 


32.43 


0.0324 


1.1 


1.15 ±0.08 


0.70 ± 0.06 


0.94 ± 0.08 


Abell 2219 


0.226 


16 


40 


19.900 


+46 


42 


41.00 


0.0178 










Abell 2390 


0.233 


21 


53 


37.115 


+17 


41 


46.41 


0.0680 


1.1 


3.90 ± 0.26 


1.86 + 0.17 


2.54 ± 0.26 


Abell 2667 


0.230 


23 


51 


39.218 


-26 


05 


03.49 


0.0165 


7.6 


2.10 + 0.12 


0.81 ± 0.08 


1.11+0.12 



The XMM-Newton observations of Abell 2219 are flared. 

Table 2. Deduced properties of all 13 galaxy clusters. Column (1): cluster name; Col. (2): electron number density of the inner 
most radial bin; Col. (3): entropy at 0.2r5oo; Col. (4): inner most radial bin of the surface brightness profile; Col. (5): cooling time 
measured at Vcen, Col. (6): cooling radius; Col. (7): rsoo; Cols. (8-10): gas mass, total mass and gas mass fraction at rsoo. 



Name 




n«o 




^cen 




tc 


eool 


'■500 








/gas,500 




10-^ 


cm ' 


keV cm^ 


kpc 




Gyr 


'•500 


Mpc 


10"Me 


10"*Me 






Abell 68 


6.4 


±0.3 


487 ± 35 


3.55 


6.6 


±0.2 


0.10 


1.21 


0.68 ± 0.07 


6.51 + 1.93 


0. 


105 + 0.060 


Abell 209 


5.7 


±0.2 


370 ± 26 


3.74 


6.4 


±0.3 


0.14 


1.15 


0.77 ± 0.08 


5.33 + 1.71 


0. 


146 ± 0.089 


Abell 267 


7.4 


±0.4 


365 ± 33 


3.29 


4.9 


±0.4 


0.15 


1.06 


0.47 ± 0.04 


4.29 ± 1.30 


0. 


109 ± 0.060 


Abell 383 


54.6 


± 1.5 


303 ± 14 


1.98 


0.5 


±0.2 


0.16 


0.98 


0.33 ± 0.04 


3.17 ±0.94 


0. 


104 ± 0.061 


Abell 773 


7.4 


±0.4 


434 ± 41 


3.15 


5.7 


±0.5 


0.14 


1.33 


1.05+0.12 


8.30 ± 2.45 


0. 


126 ± 0.070 


Abell 963 


14.3 


±0.8 


393 ± 18 


2.14 


2.6 


±0.2 


0.14 


1.14 


0.62 ± 0.07 


5.19 ± 1.52 


0. 


120 ± 0.066 


Abell 1689 


30.1 


±0.9 


505 ± 25 


1.95 


1.5 


±0.2 


0.13 


1.44 


1.05 ±0.14 


10.26 ± 2.96 


0. 


102 ± 0.060 


Abell 1763 


6.5 


±0.3 


402 ± 47 


4.00 


6.1 


±0.5 


0.13 


1.12 


0.88 ± 0.08 


4.96 ± 1.46 


0. 


178 ± 0.094 


Abell 1835 


60.8 


± 1.2 


334 ± 20 


2.50 


0.6 


±0.3 


0.19 


1.22 


1.14±0.11 


6.62 ± 1.94 


0. 


172 ± 0.092 


Abell 2218 


5.8 


±0.2 


444 ±51 


3.26 


8.0 


±0.5 


0.12 


1.07 


0.62 ± 0.06 


4.18 ± 1.27 


0. 


147 ± 0.085 


Abell 2390 


40.1 


± 1.6 


486 ± 34 


2.31 


0.9 


±0.3 


0.14 


1.29 


1.21 ±0.16 


7.67 ± 2.28 


0. 


158 ±0.095 


Abell 2667 


40.7 


± 1.4 


399 ± 27 


2.33 


0.8 


±0.3 


0.15 


1.19 


0.77 ± 0.09 


6.02 ± 1.74 


0. 


128 ± 0.073 


Mean 
























0.13 + 0.01 
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Table 3. Generalized NFW model fit. Column (1): cluster name; Cols. (2,3): characteristic density and scale of the halo of the 
generalized NFW fit; Col. (4): slope parameter of the generalized NFW fit. 



Name 


Or 




' s 












lO'^Mo Mpc"^ 




Mpc 










Abell 68 


82 + 3 


0. 


203 + 005 


0.000 


+ 


0, 


,100 


AhHl 900 


38 + 5 


0. 


989 + 018 


0.000 




0, 


.063 


Abell 267 


82 ± 5 


0. 


184 + 0.008 


0.000 


+ 


0. 


,100 


Abell 383 


108 ± 15 


0. 


1 27 + 008 


0.518 


+ 


0, 


044 


Abell 773 


53+2 


0. 


262 + 0.006 


0.000 


+ 


0. 


,100 


Abell 963 


173 ± 33 


0. 


126 + 0.010 


0.000 


+ 


0. 


,090 


Abell 1689 


64± 14 


0. 


220 ± 0.021 


0.608 


+ 


0, 


,070 


Abell 1763 


69 ±6 


0. 


,187 + 0.010 


0.000 


± 


1, 


,000 


Abell 1835 


275 ±11 


0. 


,117 + 0.002 


0.000 


+ 


0, 


.100 


Abell 2218 


123 + 4 


0. 


149 + 0.003 


0.000 


+ 


1, 


.000 


Abell 2390 


128 + 20 


0. 


,168 + 0.011 


0.324 


± 


0, 


,055 


Abell 2667 


165 ± 14 


0. 


,142 + 0.005 


0.037 


+ 


0, 


,035 



Table 4. Cluster global temperatures and metalficities. Column (1): cluster name; Col. (2,3): volume averaged radial temperature 
profile of 0.1-0.5 rsoo and 0.2-0.5 rsoo, respectively; Col. (4-9): spectral measured temperatures and metaUicities of the annuli of 
R < 2/3rt, O.lrsoo <R< 0.5r5oo and 0.2r5oo <R< 0.5r5oo, respectively. 



Name 


T'o.l-O.Srjoo 


7^0.2-0.5^500 


<2/3n 




0.2-0.5r™„ 


2<2/3r, 


2o.l-0.5r5oo 


Zo.2-0.5r5oo 




keV 


keV 


keV 


keV 


keV 


Zo 


Zo 


Zo 


Abell 68 


8.2 + 0.3 


7.7 + 0.3 


7.9 + 0.3 


7.7 + 0.3 


7.3 + 0.3 


0.17 + 0.05 


0.15 + 0.04 


0.12 + 0.05 


Abell 209 


64 + 0.2 


7.1+0.3 


6.9 + 0.2 


7.0 + 0.2 


7.1+0.3 


0.28 + 0.04 


0.23 + 0.03 


0.22 + 0.04 


Abell 267 


6.1+0.2 


6.5 + 0.4 


6.1+0.2 


6.1+0.2 


6.2 + 0.4 


0.23 + 0.05 


0.23 + 0.05 


0.25 + 0.08 


Abell 383 


4.6 + 0.1 


5.3+0.2 


4.5+0.1 


4.7 + 0.1 


4.7 + 0.2 


0.44 + 0.04 


0.29 + 0.04 


0.18 + 0.06 


Abell 773 


7.8 + 0.3 


8.1+04 


7.9 + 0.2 


7.8 + 0.3 


8.3+0.4 


0.30 + 0.03 


0.26 + 0.04 


0.30 + 0.06 


Abell 963 


6.5+0.2 


6.3+0.2 


6.4 + 0.2 


6.5 + 0.2 


6.5 + 0.2 


0.31+0.04 


0.30 + 0.03 


0.28 + 0.05 


Abell 1689 


9.0 + 0.2 


8.4 + 0.2 


9.0 + 0.1 


9.0 + 0.2 


8.5 + 0.2 


0.26 + 0.02 


0.24 + 0.02 


0.26 + 0.04 


Abell 1763 


7.1+0.2 


6,3 + 0.3 


6.3 + 0.2 


6.1+0.2 


5.8 + 0.3 


0.27 + 0.05 


0.23 + 0.04 


0.20 + 0.05 


Abell 1835 


7.6 + 0.2 


8.0 + 0.3 


7.2 + 0.1 


8.1+0.2 


8.4 + 0.3 


0.29 + 0.02 


0.25 + 0.02 


0.23 + 0.04 


Abell 2218 


8.0 + 0.2 


74 + 0.3 


6.9 + 0.2 


6.9 + 0.2 


6.6 + 0.3 


0.25 + 0.04 


0.18 + 0.21 


0.21 + 0.04 


Abcll 2390 


9.7 ± 0.4 


10.6 + 0.6 


10.1 +0.3 


10.6 + 0.4 


11.6 + 0.6 


0.35 ± 0.04 


0.28 ± 0.04 


0.24 ± 0.06 


Abcll 2667 


7.0 ± 0.2 


7.3 ± 0.3 


6.6 ± 0. 1 


7.2 ± 0.2 


7.0 ± 0.3 


0.34 ± 0.03 


0.30 ± 0.03 


0.29 ± 0.04 
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Table 5. Power law, Y = Yq X^, parameterized X-ray scaling relations. 



X 


Y 




Yn 
^ 




Scatter 


SI amnl p 












X 


Y 




keV 


^o-^-iw pc-14/3 
keVcm2 


te) 


J0208±0.26 

J02-00+0.16 
lQl.96±0.22 
1q2.04±0.18 


0.90 ± 0.30 

0.81 ±0.19 
SQ + 9'5 

\J.O-7 IE \J.Z.^ 

0.71 + 0.21 


z 


I 


REFLEX-DXL 
pilot LoCuSS sample 

iwjii V-'V-'V-' oiiuaaiiiuiC' 

CCC subsample 


r 
'•500 




y^cfij \ WkeV ) 


J03-52±0.03 

J 03.50+0.02 
103.49±0.02 
103-52±0.02 


1.00 ±0.07 

1.01 ±0.04 
97 + 05 
1.10 + 0.05 


— 


— 


REFLEX-DXL 

pilot LoCuSS sample 
non-OCO siihsainnlp 

CCC subsample 


^0-2-0.5r5oo 

keV 




r 


10l3.69±0.03;j-^l 


1.5 (fixed) 


0.05 


0.15 


pilot LoCuSS sample 








1013.72*0.04;^-^ 
10l3.91±0.38;,^-^l 
10l3.41±0,53;j-^l 


1.5 (fixed) 
1.24 ±0.44 
1.87 + 0.63 


0.05 
0.06 
0.07 


0.20 
0.15 
0.19 


non-CCC subsample 
pilot LoCuSS sample 
non-CCC subsample 



keV Mo \2 




10l2.54±0.04;j-^l 


1.8 (fixed) 


0.06 


0.10 


pilot LoCuSS sample 






inl2.54+0.05l,-l 


1.8 (fixed) 


0.06 


0.11 


non-CCC subsample 






10H-72±0.70;,-^l 


2.76 ±0.81 


0.05 


0.13 


pilot LoCuSS sample 






10n.65±0.98;j-^l 


2.86+1.17 


0.05 


0.13 


non-CCC subsample 



7'o.2-0.5.5oo 

keV 




rcorr / . \-0.5 
ergs-1 '^^^> \A,,oj 


iri42.3l±0.04L-l 
iU 

1042.30±0.05;,-^1 
1041.87±0.71;,-^1 
1042.17±0,86;j-^l 


2.60 (fixed) 

2.60 (fixed) 
3.11 ±0.82 
2.76+1.01 


0.05 

0.05 
0.04 
0.05 


0.12 

0.13 
0.13 
0.14 


pilot LoCuSS sample 

non-CCC subsample 
pilot LoCuSS sample 
non-CCC subsample 


^0-2-0.5^500 

keV 




ergs-i \a,,o) 


1042.39±0.04;,-1 

1042.38*0.05;,-^ 
1041.84±0.70;j-^l 
1041.95±0.90;j-^l 


2.98 (fixed) 

2.98 (fixed) 
3.62 + 0.80 
3.49 + 1.06 


0.04 

0.05 
0.04 
0.04 


0.12 

0.13 
0.13 
0.14 


pilot LoCuSS sample 

non-CCC subssimple 
pilot LoCuSS sample 
non-CCC subsample 






/corr / . \-0.5 
0.1-2.4keV c/•_^-l I 'ic.z 1 
ergs-1 ^A,,oj 


1018.65+0.04;,-^! 

1018.61+0.04;,-^! 
10!8.!2±4.74;,-^! 

1022.85+351;,-^! 


1.73 (fixed) 

1.73 (fixed) 
1.77 + 0.32 
1.45 + 0.23 


0.14 

0.19 
0.14 
0.19 


0.11 

0.12 
0.11 
0.12 


pilot LoCuSS sample 

non-CCC subsample 
pilot LoCuSS sample 
non-CCC subsample 






ergs-1 (aI,o) 


1015.16+0.04 

1015.11+0.04 

1014.16+5.39;,-^! 

1019.88+3.19;,-^! 


1.99 (fixed) 

1.99 (fixed) 
2.06 ± 0.36 
1.67 ±0.21 


0.14 

0.19 
0.14 
0.19 


0.12 

0.13 
0.12 
0.13 


pilot LoCuSS sample 

non-CCC subsample 
pilot LoCuSS sample 
non-CCC subsample 
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r (Mpc) r (Mpc) 



Fig. 1. De-projected temperature profiles (upper panels), surface brightness profiles (middle panels) and mass profiles (lower 
panels). The temperature profiles are approximated by the parameterization T(r) - Tj, exp[-(r-Ti)^/T2] + T(,{1 + r^/T^)~^^ + Tj 
crossing all the data points (solid). The vertical line denotes 0. 1 rjoo- The surface brightness profiles are parameterized by double-jS 
model. The observed mass profiles are fitted by generalized NFW model. 




Fig. 2. Left: Scaled radial temperature profiles for the pilot LoCuSS sample. The CCCs are in triangles. See the electronic edition 
of the Journal for a color version. Right: An average temperature profile of the pilot LoCuSS sample (black, solid) compared to 
the temperature profile ranges in Markevitch et al. (1998, grey, solid), Vikhlinin et al. (2005, grey, dashed), Zhang et al. (2006, 
black, dashed) and Pratt et al. (2007, black, dash-dotted). 
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Fig. 3. Scaled metallicity profiles. The values in the last bins for some clusters are fixed to 0.3 Z© in the spectral fitting which 
have no error bars. The CCCs are in triangles. The colors have the same meaning as those in Fig. |2] Open and filled boxes denote 
the average of the scaled metallicity profiles of the CCCs and non-CCCs, respectively, in De Grandi & Molendi (2002, grey). 




am 0.1 



Fig. 4. Surface brightness profile fits scaled according to the empirical scaling, Sx °^ The colors have the same meaning as 
those in Fig.|2l 




Fig. 5. Cooling time profiles. The CCCs are in triangles. The horizontal line denotes the age of the Universe at z ~ 0.2. The colors 
have the same meaning as those in Fig. |2] 
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Fig. 6. Entropy at 0. 1 r2oo (left) and 0.2r5oo (right) vs. temperature for the pilot LoCuSS sample (crosses) in which the CCCs are in 
triangles. Nearby clusters in Ponman et al. (2003, boxes) and the REFLEX-DXL sample (diamonds) are shown for comparison. 
The lines denote the best fits for the pilot LoCuSS sample (solid) and the REFLEX-DXL sample in Zhang et al. (2006, dashed). 
The colors have the same meaning as those in Fig.|2l 




To.2-o.5r„. (keV) 



Fig. 7. Entropy at O.Srsoo (left) and 0.3r2oo (right) vs. temperature for the pilot LoCuSS sample (crosses), the REFLEX-DXL 
sample (diamonds) and the sample in Pratt et al. (2006, boxes). The Une denote the best fit in Pratt et al. (2006, including only 
the E{z) correction). The pilot LoCuSS CCCs are in triangles. The colors have the same meaning as those in Fig.|2] 




Fig. 8. Left: Scaled entropy profiles for the pilot LoCuSS sample and the combined best fit. The shadow denotes the range of the 
entropy profiles of the nearby clusters in the same temperature range in Ponmen et al. (2003). Right: Scaled entropy profile fits 
for the pilot LoCuSS sample compared to the REFLEX-DXL sample in Zhang et al. (2006, grey, thin) and the sample in Pratt et 
al. (2006, black, hatched). The CCCs are in triangles. Tio denotes To 2-0.5^500/ 10 keV. The colors have the same meaning as those 
inFig.|2] 
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Fig. 9. Scaled total mass profiles (left) and gas mass profiles (right). The CCCs are in triangles. The colors have the same meaning 
as those in Fig.|2] The error bars (a few per cent) of the gas mass profiles are too small to be seen in this plot. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



t'- bol bol)/^ bol 



Fig. 10. Normalized cumulative cluster number count as a function of the fraction of the total bolometric luminosity attributed 
by the < 0.2r5oo region for the pilot LoCuSS sample at z ~ 0.2 compared to the REFLEX-DXL sample at z ~ 0.3 in Zhang et al. 
(2006). 




Fig. 11. Mass temperature relation (left) and gas mass temperature relation (right). The CCCs are in triangles. The colors have 
the same meaning as those in Fig. |2] The best fit power law of the pilot LoCuSS sample and REFLEX-DXL sample used the 
fixed slope of 1.5 for the M-T relation and 1.8 for the Mg^^-T relation. 
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Fig. 12. X-ray luminosity in the 0.1-2.4 keV band vs. temperature (left) and bolometric luminosity vs. temperature (right). The 
CCCs are in triangles. The clusters in Kotov & ViMilinin (2005) are in diamonds. The colors have the same meaning as those 
in Fig. |2] The best fit power law of the pilot LoCuSS sample and REFLEX-DXL sample used the fixed slope of 2.6 for the 
4kev-^ relation and 2.98 for the Ll°J^-T relation. 




Fig. 13. X-ray luminosity in the 0.1-2.4 keV band vs. total mass (left) and bolometric luminosity vs. cluster total mass (right). 
The CCCs are in triangles. The colors have the same meaning as those in Fig. |2l The best fit power law of the pilot LoCuSS 
sample used the fixed slope of 1.73 for the i'o™l2 4keV~^ relation and 1.99 for the L^™-M relation. The fixed slope of 1.3 was 
used for the REFLEX-DXL sample in Zhang et al. (2006). 




Projected X-ray mass (lO^^Mg) 



Fig. 14. Strong lensing mass vs. projected X-ray mass at rasoo for the 805 subsample (boxes) and weak lensing mass vs. projected 
X-ray mass at r™yy for the B06 subsample (circles), respectively. The symbols for the CCCs are solid. The colors have the same 
meaning as those in Fig.|2] 
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Lensing to X-ray projected mass ratio 



Fig. 15. Left: r25oo vs. strong lensing to X-ray projected mass ratio for the 805 subsample (boxes) and r™QQ vs. weak lensing to 
X-ray projected mass ratio for the B06 subsample (circles), respectively. Right: substructure fraction (defined in Smith et al. 
2005) vs. strong lensing to X-ray projected mass ratio at r25oo for the S05 subsample (boxes) and weak lensing to X-ray projected 
mass ratio at r^QQ for the B06 subsample (circles), respectively. The symbols for the CCCs are solid. The colors have the same 
meaning as those in Fig.|2] 



Projected gas mass fraction using lensing masses 



Fig. 16. Projected X-ray gas mass fraction vs. projected gas mass fraction using strong lensing (boxes) and weak lensing (circles) 
masses, for the S05 subsample and B06 subsample, respectively. The symbols for the CCCs are solid. The colors have the same 
meaning as those in Fig. |2] 




Fig. 17. Flat fielded XMM-Newton MOSl flux images in the 0.7-2 keV band superposed with the CFH12k optical luminosity 
weighted galaxy density contours. 
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Fig. 18. X-ray luminosity in the 0.1-2.4 keV band vs. strong lensing mass within r25()o for the 805 subsample (boxes), and weak 
lensing mass within rjgo for the B06 subsample (circles), respectively. The symbols for the CCCs are solid. The colors have the 
same meaning as those in Fig.|2l The slope of the fitting is set to 1.73. The grey lines denote the best fit power law for the 805 
subsample (dashed) and B06 subsample (solid). The black line denotes the X-ray luminosity-mass relation for the pilot LoCu88 
sample at z ~ 0.2. 
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Table A.l. Observational information of the XMM-Newton data. 



Name 


Date 


Id 


Filter 


Frame 


Net exposure (ks) 










MOS 


r>n 


MOSl 


M0S2 


nn 


Abell 68 


2006-03-14 


0084230201 


Medium 


FF 


EFF 


24.94 


23.81 


18.15 


Abell 209 


2001-01-15 


0084230301 


Medium 


FF 


EFF 


17.33 


16.09 


12.80 


Abell 267 


2002-01-02 


0084230401 


Medium 


FF 


EFF 


16.06 


15.74 


10.96 


Abell 383 


2002-08-17 


0084230501 


Medium 


FF 


EFF 


28.07 


28.05 


21.50 


Abell 773 


2001-04-26 


0084230601 


Medivmi 


FF 


EFF 


13.01 


14.68 


15.90 


Abell 963 


2001-11-02 


0084230701 


Medium 


FF 


EFF 


23.82 


24.69 


17.92 


Abell 1689 


2001-12-24 


0093030101 


Thin 


FF 


EFF 


33.67 


33.34 


29.24 


Abell 1763 


2002-12-13 


0084230901 


Medium 


FF 


EFF 


12.32 


12.01 


9.34 


Abell 1835 


2000-06-28 


0098010101 


Thin 


FF 


FF 




25.27 


24.73 


Abell 2218 


2002-09-28 


0112980101 


Thin 


FF 


EFF 


16.69 


16.91 


13.79 


Abell 2219 


2002-06-14 


0112231801 


Medium 


FF 


EFF 










2002-06-24 


0112231901 


Medium 


FF 


EFF 








Abell 2390 


2001-06-19 


0111270101 


Thin 


FF 


FF 


10.27 


10.03 


8.85 


Abell 2667 


2003-06-21 


0148990101 


Medium 


FF 


FF 


22.21 


22.87 


14.10 



The MOSl data of Abell 1835 are in window mode, which cannot be used for this work. The XMM-Newton observations of Abell 2219 are 
flared. 
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Fig.B.l. See the caption in Fig.[T] 
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Fig.B.2. See the caption in Fig.[T] 
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Fig.B.3. See the caption in Fig.[T] 
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Fig.B.4. See the caption in Fig.[T] 
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Fig. B.5. See the caption in Fig.[T] 
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Fig.B.6. See the caption in Fig.[n] Abell 773 has no CFH12k data. 




Fig.B.7. See the caption in Fig.[ni Abell 2667 has no CFH12k data. 



